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SECTION  I 
INTRODUCTION 

This  final  report  summarizes  the  results  which  have  been  obtained 
under  ARO  Contract  DAAG  29-77-C-00I0.  The  purpose  of  this  research  has 
been  the  investigation  of  the  infrared  water  vapor  continuum.  Although 
a  great  deal  of  data  existed  prior  to  this  research  effort,  a  careful 
analysis  of  the  continuum  absorption  data  had  not  been  performed.  Empir¬ 
ical  interpretations  of  experimental  data  provided  convenient  recipes 
for  calculating  continuum  absorption  values  in  specific  spectral  regions 
and  under  certain  conditions  of  pressure,  water  vapor  content,  and  tem¬ 
perature.  However,  predictions  of  continuum  absorption  in  uninvestigat¬ 
ed  spectral  regions  or  for  physical  conditions  outside  the  body  of  exist¬ 
ing  experimental  data  were  not  possible. 

Data  which  have  been  collected  under  this  research  program  include 
room  temperature  and  elevated  temperature  transmittance  values  of  water 
vapor  samples  broadened  with  nitrogen  and  with  air  (80%  N^,  20%  0^).  A 
line  tunable,  frequency-stabilized  CO^  laser  was  used  as  the  radiation 
source  for  long-path  absorption  cell  measurements  and  for  nonresonant 
spectrophone  measurements.  Also,  high-resolution  Fourier  transform 
spectroscopy  was  employed  to  record  pure  and  pressure-broadened  water 
vapor  spectra  under  a  variety  of  conditions.  These  spectra  are  discuss¬ 
ed  later  in  this  report  and  will  ultimately  be  made  available  as  an  at¬ 
las  of  H^O  absorption  lines. 

We  believe  that  the  infrared  water  vapor  continuum  is  caused  en¬ 
tirely  by  the  wings  of  strong  absorption  lines  of  the  H^O  molecule.  The 
strongest  absorption  lines  of  the  water  molecule  are  located  between 
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200  cm-*  and  400  cm~*.  The  wings  of  these  absorption  lines  contribute 
significantly  to  the  continuum  absorption  even  in  the  8-14  pm  atmosphere 
transmittance  window,  but  have  little  effect  in  the  3. 5-4. 2  pm  spectral 
region. 

A  detailed  theory  of  the  spectral  line  shape  for  pure  and  pressure- 
broadened  H^O  molecules  has  been  formulated.  The  resulting  non-Lorent- 
zian  line  shape  predicts  the  negative  temperature  dependence  observed  in 
the  20  pm  region. 

Section  II  of  this  report  contains  a  review  of  the  early  work  con¬ 
cerning  water  vapor  continuum  absorption.  Section  III  summarizes  the 
theoretical  developments  which  led  to  the  formulation  of  the  spectral 
line  shape  for  H^O  molecules.  Section  IV  discusses  the  Fourier  trans¬ 
form  spectra  which  were  collected.  These  spectra  should  prove  useful  to 
those  interested  in  locating  the  weak  H^O  absorption  lines  within  the 
atmospheric  transmittance  windows. 
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SECTION  II 

THE  WATER  VAPOR  CONTINUUM 

The  transmission  of  infrared  radiation  through  the  atmosphere  con¬ 
tinues  to  be  a  subject  of  considerable  interest.  Knowledge  of  the  vari¬ 
ous  attenuation  processes  due  to  both  molecular  absorption  and  aerosol 
scattering  is  important  to  many  different  scientific  disciplines.  In 
meteorologic  studies  concerning  the  energy  budget  for  the  surface  of  the 
earth,  for  example,  contributions  from  radiation  transmitted  through  the 
atmosphere  to  the  surface,  as  well  as  energy  radiated  away  from  the  sur¬ 
face,  must  be  known.  The  transmission  functions  which  describe  these 
processes  are  dependent  on  molecular  constituents,  particularly  H^O  and 
CO^,  and  also  on  the  aerosol  concentrations  in  the  atmosphere. 

Infrared  transmission  studies  are  also  needed  to  determine  optimal 
spectral  regions  for  infrared  imaging  systems  and  the  evaluation  of  in¬ 
frared  sensor  performance  under  various  atmospheric  conditions.  Imaging 
and  sensor  systems  must  prove  reliable  in  the  field  under  all  possible 
circumstances. 

Another  area  in  which  precise  information  on  the  attenuation  mech¬ 
anisms  is  needed  is  the  study  of  the  propagation  of  infrared  laser  radi¬ 
ation.  The  transmission  of  laser  radiation  through  the  atmosphere  has 
many  applications.  When  considering  long  path  lengths  through  the  atmos¬ 
phere  for  projects  such  as  communications,  laser  ranging,  or  energy 
transfer,  even  small  attenuations  per  meter  can  become  serious. 

Early  workers  in  infrared  transmission  studies  recorded  the  spec¬ 
tral  transmittance  of  the  atmosphere  to  infrared  solar  radiation.  S.P. 
Langley  LI, 2]  published  several  infrared  solar  spectra.  From  his  work. 
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he  concluded  that  there  was  practically  complete  transmission  at  10.7cm. 

A.  Adel  L3,4]  repeated  much  of  Langley's  work  and  extended  the  solar 
spectrum  to  24  pm.  While  at  The  Ohio  State  University,  M.V.  Migeotte 
[5,6]  investigated  the  infrared  solar  spectrum  from  2  pnr.  to  12  pm.  This 
work  was  continued  in  the  early  1950 ' s  by  J.H.  Shaw  and  co-workers  [7]. 

The  emphasis  of  this  early  work  was  on  the  identification  of  atmos¬ 
pheric  gases  and  the  accurate  determination  of  absorption  features. 

Many  high-resolution  maps  of  the  solar  spectrum  were  published  which  col¬ 
lectively  cover  the  region  from  1  pm  to  350  pm. 

In  1938,  W.M.  Elsasser  L8]  suggested  the  existence  of  a  water  vap¬ 
or  continuum  absorption  in  the  8  pm  to  12  pm  region  which  he  attributed 
to  extreme  wings  of  very  strong  absorption  lines  outside  this  region. 

Yates  and  Taylor  [9]  made  measurements  over  horizontal  paths  using  low 
resolution  (10  cm*1)  spectroscopy  and  found  greater  attenuation  than 
expected  from  selective  absorption  alone.  The  existence  of  continuum 
absorption  near  10  pm  has  also  been  confirmed  in  solar  spectra  by  A.  Adel 
[10],  R.  Anthony  [11],  and  W.T.  Roach  and  R.M.  Goody  [12].  In  the  early 
1960's,  K.J.  Bignell  and  co-workers  L 1 3 ]  studied  continuum  absorption  in 
the  10  pm  region  at  1  to  2  cm”1  resolution  using  the  sun  as  a  source, 
and  in  the  20  pm  region  using  horizontal  atmospheric  paths  and  a  Nernst 
source.  They  guessed  that  the  temperature  dependence  of  the  continuum 
absorption  would  be  negligibly  small.  Therefore,  their  data  are  report¬ 
ed  without  indication  of  the  temperature.  In  a  later  paper,  however, 
Bignell  [14]  reported  a  strong  negative  dependence  of  the  continuum  absorp¬ 
tion  coefficient  on  the  temperature,  and  also  reported  contributions  to 
this  absorption  coefficient  from  a  self-broadening  term  and  a  foreign 
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broadening  term. 

D.E.  Burch  [15]  has  made  laboratory  studies  of  the  water  vapor  con¬ 
tinuum  in  the  8  pm  to  12  pm  region  and  also  in  the  4  pm  region  [16].  In 
the  8  pm  to  12  pm  region  he,  too,  found  that  the  continuum  absorption  co¬ 
efficient  decreased  rapidly  with  temperature.  Furthermore,  he  found  that 
the  calculated  values  for  the  continuum  absorption  coefficient  based  on 
the  Lorentz  line  shape  for  the  strong  water  vapor  lines  were  lower  than 
the  experimentally  measured  absorption  coefficients  at  all  wavelengths  in 
the  8  pm  to  12  pm  region. 

P.  Varanasi,  S.  Chou,  and  S.S.  Penner  [17]  have  suggested  that  the 
continuum  absorption  in  this  region  is  due  to  absorption  by  water  dimer 
molecules.  The  absorption  coefficient  from  this  mechanism  should  in¬ 
crease  with  pressure  in  the  same  way  as  absorption  due  to  the  extreme 
wings  of  strong  lines.  The  basis  for  the  dimer  hypothesis  is  the  strong 
temperature  dependence  which  is  exhibited  by  the  continuum  absorption. 

CO^  laser  measurements  of  the  8  pm  to  12  pm  water  vapor  continuum 
have  been  made  at  this  laboratory  by  J.H.  McCoy,  D.B.  Rensch,  and  R.K. 
Long  [18].  Recent  measurements  have  been  made  by  F.S.  Mills,  R.K.  Long, 
and  E.K.  Damon  [19].  These  measurements  involved  the  use  of  a  multiple- 
traversal  absorption  cell  to  achieve  path  lengths  of  1  km.  Samples  of 
pure  water  vapor,  as  well  as  water  vapor-air  mixtures,  were  studied  at 
several  output  frequencies  of  the  CO.,  laser.  At  \  =  10.59  pm,  McCoy 
found  that  the  self-broadening  contribution  to  the  continuum  absorption 
coefficient  was  about  200  times  stronger  than  the  contribution  due  to 
foreign  gas  broadening  [18]. 

Studies  of  the  absorption  of  C0^  laser  radiation  in  the  atmosphere 
are  also  being  performed  in  the  Soviet  Union.  Absorption  coefficients 


for  the  continuum  absorption  at  the  centers  of  the  P(1 2)  to  P(2 6)  laser 
lines  in  the  10.4  pm  band  of  CO2  have  been  studied  by  I. I.  Ippolitov  [20] 
and  by  V.N.  Aref'ev,  et  al . ,  [21].  A  complete  review  of  the  influence 


of  molecular  absorption  on  the  propagation  of  CO^  laser  radiation  in  the 
atmosphere  is  given  by  T.G.  Adiks,  V.N.  Aref'ev,  and  V.I.  Dianov-Klokov 
[22]. 

Laboratory  measurements  of  the  infrared  continuum  can  be  divided 
into  two  groups,  depending  on  whether  a  laser  or  blackbody  radiation 
source  is  used.  The  advantage  of  the  laser  measurement  is  that  more  ac¬ 
curate  transmittance  values  are  obtainable  due  to  the  high  power,  and 
excellent  frequency  and  amplitude  stability.  The  major  disadvantage  is 
that  the  radiation  is  available  at  only  a  few  fixed  frequencies.  Inform¬ 
ation  on  the  continuum  absorption  at  a  particular  frequency  could  be  ob¬ 
scured  by  a  dominating  local  absorption  line.  Local  and  continuum  ab¬ 
sorption  have  different  partial  pressure  dependencies,  so  that  this  may 
be  used  as  a  test  for  local  effects. 

Broadband  spectroscopic  laboratory  measurements  depend  on  comparing 
two  transmittance  spectra.  The  first  spectrum  is  recorded  with  the  ab¬ 
sorption  cell  evacuated.  The  second  spectrum  is  recorded  with  the  water 
vapor  sample  in  the  cell,  and  determination  of  the  100%  transmittance 
curve  is  sometimes  difficult. 

Pressure  Dependence 

The  8  pm  to  12  pm  atmospheric  transmission  window  plays  an  import¬ 
ant  role  in  the  energy  budget  of  the  surface  of  the  earth,  and  in  the 
heat  balance  of  the  atmosphere  since  the  spectral  radiance  of  a  blackbody 
near  296K  is  a  maximum  in  this  region.  Infrared  imaging  and  sensor 
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systems  also  make  extensive  use  of  this  spectral  region.  Attenuation  of 
radiation  by  continuum  absorption  in  this  region,  therefore,  is  import¬ 
ant. 

A  comprehensi ve  study  of  the  water  vapor  continuum  will  provide 
the  necessary  data  for  systems  design  and  should  also  provide  an  impetus 
for  further  theoretical  understanding  of  the  mechanism  of  continuum  ab¬ 
sorption. 

A  more  recent  report  on  the  water  vapor  continuum  in  the  8  pm  to 
12  pm  region  is  by  R.E.  Roberts,  L.M.  Biberman,  and  J.E.A.  Selby  [ 233 . 
They  have  accumulated  much  of  the  previously  published  data,  including 
laser  measurements  from  this  laboratory,  and  also  unpublished  spectro¬ 
scopic  data  by  D.E.  Burch,  in  an  attempt  to  evaluate  the  accuracy  of  the 
LOWTRAN  continuum  model  of  J.E.A.  Selby  and  R.A.  McClatchey  [24] . 

The  amount  of  radiant  energy  at  frequency  v  transmitted  through 
an  absorber  of  length  L  (in  cm)  can  be  expressed  as 

T(v)  =  exp(-k(v)L)  (1) 

where  k(v)  is  the  absorption  coefficient  (in  cm’^)  at  frequency  v.  This 
absorption  coefficient  is  a  combination  of  the  attenuation  due  to  molecu¬ 
lar  absorption  and  radiation  loss  due  to  aerosol  scattering.  In  this 
discussion,  however,  it  is  assumed  that  the  contribution  from  aerosols 
is  small.  From  Equation  (1),  the  absorption  coefficient  can  be  written 

k(v)  =  — j~  In  T (v )  .  (2) 

Before  this  topic  is  developed  in  great  detail,  it  is  important 
to  compare  this  starting  point  with  methods  used  by  other  workers.  Burch 
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[153  used  a  different  starting  point  by  writing 
T(v)  =  exp(-k'(v)u) 


(3) 


as  his  first  equation.  Here  k'(v)  is  the  absorption  coefficient  in 

O  1 

(molecule  •  cm  )  .  We  have  included  the  prime  to  distinguish  this 

coefficient  from  Equation  (2).  The  quantity  u  is  the  absorber  thickness 

_2 

express  in  molecules  •  cm  .  Burch  then  writes 


k»  -4'n  T<u>  =  ("i  ph/c"(v)Pn 


(4) 


for  the  continuum  absorption  coefficient.  This  expression  is  based 
on  the  assumption  that  the  absorption  is  due  to  the  extreme  wings  of 
collision-broadened  lines.  In  Equation  (4),  C°  (v)  and  C°  (v)  are  k * ( v ) 
from  self-broadening  and  foreign  gas  broadening  (usually  nitrogen)  mech¬ 
anisms.  Burch  demonstrated  the  linear  dependence  of  the  continuum  ab¬ 
sorption  coefficient  k'(v)  on  the  water  vapor  pressure  for  a  pure  sample 
of  H2O  (P^  =  0)  by  plotting  ^  In  T(v)  versus  q  for  data  taken  at 
v  =  844.2  cm  at  a  temperature  of  387  K.  This  plot  is  shown  in  Figure 
1. 


Although  this  linear  dependence  as  expressed  by  Burch  was  based 
upon  the  assumption  that  the  continuum  absorption  is  caused  by  extreme 
wings  of  strong  lines.  Figure  1  does  not  necessarily  substantiate  this 
mechanism.  Varanassi,  et  al . ,  [17]  have  proposed  that  continuum  absorp¬ 
tion  could  be  due  to  absorption  by  the  water  dimer  (H-^O^.  The  absorp¬ 
tion  coefficient  as  expressed  in  Equation  (4)  caused  by  absorption  by 
dimer  molecules  would  also  be  linear  in  the  water  vapor  pressure  [25]. 
This  statement  is  analyzed  in  more  detail  later  in  the  section  on  water 


MOLECULES  •  cm 


4*10 


Figure  1.  Plot  of  absorption  coefficient  k'(v)  at  v  =  844.2  cm  ‘  as  d 
fined  in  Equation  (4)  as  a  function  of  water  vapor  pressure 
for  a  sample  of  pure  water  vapor  at  T  =  387  K.  Reproduced 
from  Reference  [15]. 
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dimers.  Thus,  a  study  of  the  pressure  dependence  of  the  continuum  absorp¬ 
tion  coefficient  in  the  8  pm  to  12  pm  region  does  not  constitute  positive 
verification  of  one  absorption  mechanism  over  the  other. 

In  order  to  discuss  the  results  of  previous  research  on  the  water 
vapor  continuum,  we  will  develop  the  details  expressed  in  Equations  (1) 
and  (2).  The  absorption  coefficient  given  in  Equation  (2)  can  be  compar¬ 
ed  to  the  absorption  coefficient  in  Equation  (4)  by  noting  that 


u  ’  wH2o  '  L 


(5) 


where  wH  ^  is  the  number  density  of  water  vapor  molecules  in  units  of 
molecules  .  cm.  From  Equations  (2)  and  (4) 


^  In  T(v)  =  k(v)  =  «H  Qk'(v) 


(6) 


The  explicit  pressure  dependence  of  k(v),then,  is 

k^v)  =  wH20  [Cs(v)  PH20  +  CN^P  '  V 


(7) 


This  can  be  written 


k(v)  =  C°(v)  wH?Q 


h2o 


+  co  (P  '  PH20^ 


(8) 


In  these  two  equations,  P  is  the  total  pressure.  In  the  past,  the  assump¬ 
tion  has  been  made  that  the  ratio  cjjj(u)/C°(v)  is  independent  of  frequency. 
Roberts  points  out  that  this  might  not  be  the  case,  and  perhaps  certain 
provisions  should  be  made  to  accommodate  this  in  future  expressions.  For 
the  rest  of  this  development,  however,  it  will  be  assumed  to  be  independ¬ 
ent  of  frequency.  Thus, 
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k(v)  =  c°s(v^  wh2o  [V  +  Y  (P  -  Ph2o)  • 

Several  workers  have  measured  C°(v)  and  y  at  or  near  room  tempera¬ 
ture.  Bignell  L 143  reported  a  contribution  to  the  continuum  absorption 
coefficient  which  was  linear  in  water  vapor  pressure  that  became  equal 
to  the  foreign  broadened  term  at  a  water  pressure  of  about  15  mbar  (11.25 
torr).  From  this  information,  the  value  of  y  which  he  measured  was  near 
0.015.  Values  of  C°(v)  measured  by  Bignell  are  shown  in  Figure  5. 

Laser  absorption  studies  performed  by  our  laboratory  using  the 
P ( 20 )  laser  line  of  C02  near  10.4  pm  suggest  that  y  is  much  smaller. 
McCoy,  et  a!.,  reported  their  data  for  air  broadened  water  samples  in 
the  form  of  an  equation 


k(v  =  944.2  cm'1)  =  4.32  2  10~6  (km*  ^orr"2^  Q  P  +  193 


Comparison  between  this  equation  and  Equation  (9)  can  be  made  when 

1/y  -  1  =  193  (11-a) 

C°(944.2  cm'1)  w^Q  =  4.32  x  lO'^km^Torr'2)?^  .  (11-b) 


From  Equation  (11-a)  y  =  0.105.  This  is  a  factor  of  three  lower  than 
the  value  which  was  inferred  from  Bignell’s  work.  The  value  of 
C°(944.2  cm"1)  derived  from  McCoy's  equation  is 

C°(944.2  c.-1)  ■  .1.32_xiL-!l760,)(7.001x  l0i  296 
s  (.005)  2.687  x  10"iy  65 


-  1.85  x  10‘22cm2molecule'1atm'1 
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Burch  has  not  been  able  to  measure  y  with  any  degree  of  reliabil¬ 
ity.  Roberts  performed  a  linear  regression  to  Burch's  data  in  the  16  pm 
to  30  pm  region  and  extrapolated  his  results  to  the  8  pm  to  12  pm  region. 
Results  of  this  extrapolation  show  a  value  of  y  from  Burch's  work  to  be 
0.0008. 

Values  represented  here  for  y  differ  by  a  factor  20.  It  is  clear 
that  there  is  a  large  experimental  uncertainty  about  y.  Roberts  has 
suggested  that  until  better  measurements  are  made,  the  L0WTRAN  model 
should  eliminate  y.  Further  unknowns  about  y  include  its  frequency  de¬ 
pendence  as  discussed  earlier,  and  also  its  temperature  dependence. 

Several  studies  of  pure  water  vapor  samples  have  been  made  in  this 
laboratory.  McCoy  measured  the  transmittance  at  944.2  cm’1  and  1047.1 
cm*1  through  980  m  of  water  vapor  at  various  partial  pressures.  The  ab¬ 
sorption  coefficients  which  he  measured  at  these  frequencies  were  fitted 
to  the  equations 


k(944.2  cm'1)  =  8.39  x  10'4  P*  . 

n^U 

k ( 1047 . 1  cm'1)  =  6.54  x  10'4  P* 


(12-a) 

(12-b) 


More  recently,  F.S.  Mills  and  R.K.  Long  [19]  made  measurements  of  the 
continuum  absorption  coefficient  in  pure  water  vapor  at  the  P(20)  laser 
line  of  CO2  at  944.2  cm'1.  The  data  were  fitted  to  the  equation 


kMins(944.2  cm'1)  =  9.67  x  10'4  pg  .  (13) 

The  value  of  C°(944.2  cm"1)  obtained  from  this  fit  is  2.24  x  10"^  cm^ 
-Is  -1 

molecule  atm  A. 
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The  difference  between  the  fit  obtained  by  McCoy  and  the  fit  obtain¬ 
ed  by  Mills  is  approximately  13%.  Figure  2  shows  the  two  results. 

V.N.  Aref'ev,  et  al . ,  [21]  have  also  studied  the  transmittance  of 
pure  water  vapor  samples.  For  these  results,  a  multiple-transversal  cell 
with  a  base  path  of  50  m  was  used  to  achieve  path  lengths  of  3  km. 

The  CO.,  laser  which  was  used  as  a  source  was  not  stabilized  and  it  oscil¬ 
lated  alternately  on  the  P(16),  P(20),  P(22),  and  P(24)  transitions  in 
the  10.4  ym  band. 

Results  of  their  work  done  at  294  K  are  reproduced  in  Figure  3. 

The  dashed  curve  represents  the  results  of  McCoy,  et  al . ,  Cl8],  as  express¬ 
ed  in  Equation  (12-a).  Aref'ev  chose  to  fit  his  data  to  an  expression  of 
the  form 

T  =  exp(-(0. 031w  +  0.029w2))  (14) 

where  w  is  the  centimeters  of  precipitable  water  in  the  3  cm  path.  This 
can  be  written  for  arbitrary  path  length  in  terms  of  the  water  vapor  den- 
s i ty  a  (in  gm/m  ) 

T  =  exp(-(3.1a  +  0.87a2)  x  10‘3L)  .  (15) 

The  absorption  coefficient  in  the  form 

k  =  (3.1a  +  0.87a2)  x  10"3  (16) 

with  a  term  which  is  linear  in  the  absorber  amount  for  pure  water  samples 
is  not  in  agreement  with  the  results  of  McCoy  and  of  Mills  and  Long  in 
Figure  2.  In  the  units  chosen  by  Burch,  this  linear  term  in  Equation 
(16)  would  become  a  nonzero  intercept  to  Figure  1. 


NVER3E  KM) 
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Figure  2.  Comparison  of  the  results  of  McCoy,  et  al.,  and  the  results 
of  Mills  and  Long  for  absorption  of  the  f5 ( 20 )  laser  line  of 
C02  by  pure  water  vapor  samples. 


transmittance  as  a  function  of  precipitable  water  in  a  3  km  path.  M 
Coy’s  data  are  slightly  high  which  could  be  due  to  temperature  dif¬ 
ferences  between  the  two  experiments. 


In  light  of  this  discrepancy,  we  have  examined  the  data  of  Figure 


3  very  carefully.  Using  a  computer  graphics  technique,  we  digitized 
the  data  of  Figure  3  two  separate  times  to  create  independent  data  sets. 
We  then  ran  several  curve  fit  routines  on  these  data  sets  to  determine 


the  necessity  of  the  linear  term  in  Equation  (16). 

2 

The  two  independent  data  sets  were  fit  to  a  ^  dependence  by 

3  ^3 

converting  gm/m  to  Torr.  The  conversion  used  was  1  gm/m  =  1.033  Torr 

at  294  K.  Figures  4  (a)  and  4  (b)  show  plots  of  -(l/L)ln(T)  in  km~* 

2 

as  a  function  of  g  for  each  of  the  data  sets.  Allowance  was  made 
in  the  curve  fit  for  the  fact  that  the  best  straiqlit  line  might  not 
pass  through  the  origin  due  to  systematic  or  instrumental  errors  in  the 
experiment.  Results  of  these  fits  are  shown  in  Table  1.  Also  shown 
are  the  least  square  fits  to  McCoy's  data  and  the  data  of  Mills  and 
Long.  The  RMS  errors  of  the  two  data  sets  from  Aref'ev  indicate  a  good 
fit  to  the  data  and  we  feel  that  the  linear  term  in  Equation  (16)  is 


unnecessary. 

Comparison  of  the  values  of  C°(v)  taken  from  various  sources  is 
shown  in  Figure  5.  This  is  a  plot  of  the  self-broadening  coefficient 
as  a  function  of  frequency  in  the  8  pm  to  12  pm  region.  Measurements 
made  at  this  laboratory  on  pure  water  samples,  as  well  as  nitrogen 
broadened  water  samples,  are  shown  as  solid  points  in  the  figure.  Ori¬ 
ginal  data  published  by  Burch  in  Reference  [16]  are  high  in  the  8  pm 
region  and  show  a  great  deal  of  scatter.  Roberts  has  acquired  some 
recent  data  by  Burch  which  shows  much  less  scatter.  These  values 
are  also  shown  in  the  plot.  Bignell's  data  were  converted  from  his 

9  1 

units  of  cm  cm  to  the  units  in  Figure  5  by  multiplying  by  2.99  x 
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TABLE  1 

LEAST  SQUARE  FITS  OF  THE  ABSORPTION  COEFFICIENTS  FROM  SEVERAL  DATA 
SETS  TO  THE  FORM  k  =  a  +  b  PH  Q 


Data  Set 


01 


02 


McCoy 


Mills  and  Long 


Fi  t 


k  =  (6.19  +  1.06  P^0)  x  10 
k  -  (7.24  +  1.05  I'2  (J)  x  10 
k  =  (1.95  +  0.03  P2_Q)  x  10 


-3 


-3 


-3 


k  =  (4.94  +  0.93  P2  Q)  x  10'3 


RMS  Error 
1.3  x  10~2 

1  .3  x  l(f  2 

4.1  x  10'3 

4.0  x  10'3 


Roberts  has  carried  out  a  curve  fit  to  the  recent  data  provided 
by  Burch.  His  curve 

C°(v)  =  a  +  b*exp(-crv)  (17) 

is  not  based  on  any  theoretical  model,  but  it  does  provide  a  practical 
method  for  calculating  the  continuum  absorption  coefficient  at  any  fre¬ 
quency  v  in  the  8  pm  to  12  pm  region.  Roberts  included  Burch’s  data 
out  to  30um,  and  achieved  a  coefficient  of  determination  (r  )  of  0.992 
with  the  parameter  values  of 

a  =  1.25  x  10~33  cm3  molecule'*  atm'^ 
b  =  1.67  x  10'*9  cm3  molecule'*  atm'* 
o  =  7.87  x  103  cm 
This  fit  is  shown  in  the  plot. 
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Temperature  Dependence 

Data  presented  in  Figure  5  were  recorded  at  a  temperature  near 
296  K.  Burch  has  duplicated  these  measurements  at  several  temperatures 
above  296  K.  These  results  are  shown  in  Figure  6.  The  data  at  296  K 
shown  here  are  the  original  data  which,  according  to  Figure  5,  are 
higher  than  his  recent  data  in  the  8  Mm  region.  From  Figure  6,  the 
temperature  dependence  of  C°(v)  at  several  frequencies  can  be  found. 

We  have  examined  these  curves  at  850  cm"*,  900  cm-*,  and  1000  cm-*. 

The  three  data  points  at  each  frequency  were  then  fit  to  an  equation 
of  the  form 

C°(v)  =  a  exp(T0/T)  .  (18) 

An  equation  of  this  form  is  discussed  briefly  by  Roberts,  and  it  has 
interesting  interpretations  with  respect  to  the  dimer  hypothesis. 

At  850  cm-'*’,  no  correction  of  the  296  K  data  was  made.  At  900 
cm-*  and  1000  cm--*,  two  data  sets  were  analyzed  at  each  frequency. 

The  first  set  included  the  uncorrected  data  at  296  K,  and  the  second 
set  included  corrected  data  according  to  the  curve  fit  in  Figure  5.  No 
change  of  the  elevated  temperature  data  was  made.  The  temperature  pro¬ 
files  for  C°(v}T)  at  these  frequencies  are  plotted  in  Figure  7. 

Results  of  the  curve  fit  of  Equation  (18)  are  shown  in  Table  2. 

p 

The  high  values  for  the  coefficient  of  determination  (r  )  are  due 
to  the  fact  that  only  three  data  points  have  been  used  in  each  fit. 

The  value  of  a  is  not  expected  to  be  consistent  at  these  different  fre¬ 
quencies  since  it  contains  the  frequency  dependence  of  C°(v,T).  To  a 
first  approximation,  however,  the  value  of  T  should  be  independent  of 
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WAVENUMBER  (  cm  ) 

Figure  6.  The  self-broadening  coefficient  C°(v)  at  several  temperatures.  The  dashed  line  at  296  K 
represents  the  curve  fit  to  recent  data  in  Figure  6. 
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frequency  in  the  entire  8  pm  to  12  pm  region. 

If  the  old  data  in  Table  2  are  used,  the  average  value  of  T0  is 
found  to  be  1718  K.  The  TQ  values  from  the  more  recent  data  by  Burch 
seem  to  be  much  more  consistent  with  the  1Q  value  at  850  cm,  where  both 
old  and  new  data  are  in  agreement.  The  average  value  of  TQ  derived 
from  the  new  data  is  1638  K.  Using  this  value,  an  equation  which  is 
normalized  to  room  temperature  data  can  be  written  as: 

C°0\T)  =  Cs(v)  exp  1638  K  4 ' *  (19) 

C°(v)  is  the  self-broadening  coefficient  at  296  K. 

V.N.  Aref'ev  carried  out  measurements  of  the  continuum  absorption 
coefficients  at  the  centers  of  the  P(16),  P ( 20 ) ,  P ( 22 ) ,  and  P ( 24 ) 

C0?  laser  lines  in  the  10.4  pm  band.  As  mentioned  before,  their  laser 
oscillated  alternately  on  each  of  the  lines,  and  no  special  steps  were 
taken  to  separate  the  results  obtained.  Measurements  were  carried  out 
at  temperatures  of  293  K,  323  K,  and  353  K,  as  reported  by  Adiks,  et  al . , 
[22]  on  pure  water  samples.  Results  of  this  work  were  reported  in  the 
form  of  an  equation 


k  =  1.76  x  10'3a  +  0.42  x  10‘6a2  exp 


(20) 


In  this  equation,  a  is  the  amount  of  water  vapor  in  units  of  gm  •  m~ 
and  k  is  in  km"3.  Since  frequency  components  were  not  separated,  this 
expression  is  an  average  over  the  four  laser  lines.  In  comparing  this 
work  with  other  laser  measurements,  it  is  assumed  that  Equation  (20)  is 
valid  at  the  P(20)  laser  lines  position.  The  value  of  -AH/R  reported 


was  2273  K. 
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The  presence  of  a  term  which  is  linear  in  the  absorber  amount  has 
been  questioned  earlier  in  this  report.  The  fact  that  it  is  now  assumed 
to  be  independent  of  temperature  makes  this  term  even  more  suspect. 

An  equation  of  the  form  shown  in  Equation  (20)  is  not  in  agreement  with 
Equations  (8),  (10),  and  (19).  Furthermore,  it  does  not  even  agree  well 
with  Equation  (16),  as  can  be  seen  by  using  T  =  294  K  in  Equation  (20) 
to  write 

k  =  (1.76a  +  0.96a2)  x  10'3  .  (21) 

Other  workers  have  estimated  the  size  of  the  exponential  factor  TQ 
in  Equation  (18).  Roberts  has  summarized  these  results  in  his  report  on 
the  water  vapor  continuum  [23].  Rather  than  reporting  the  value  which 
Aref'ev  and  Dianov-Klolov  found  as  2273  K,  Roberts  fit  their  data  to  an 
equation  of  the  form  given  in  Equation  (16).  Results  of  this  fit  give 
Tq  =  1810  K.  Table  3  summarized  values  of  the  exponential  constant  found 
by  several  workers.  In  all  of  this  work,  errors  of ± 200  K  are  not  uncom¬ 
mon. 

TABLE  3 

COMPARISON  OF  VALUES  OF  T  IN  EQUATION  (18)  DETERMINED  FROM 

SEVERAL  WORKERS 


Reference 

T0  (K) 

Burch  [15]  (old  data) 

1718 

Burch  (new  data) 

1638 

Aref'ev  [22] 

2273 

Aref'ev* 

1810 

Bignell  [14] 

1800 

Varanasi,  et  al.,  [17] 

2000 

♦Using  Ay?f'ev  data  and  fitting  to  Equation  (19) 
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A  Case  for  Water  Dimer  Absorption 

Burch  found  that  the  water  vapor  continuum  absorption  coefficients 
in  the  8  urn  to  12  urn  region  are  significantly  larger  than  the  calculated 
values  based  on  extreme  wings  of  strong  Lorentz  water  lines.  Figure 
8  shows  that  the  discrepancy  is  a  full  order  of  magnitude.  The  simple 
Lorentz  line  shape  has  the  form 


k(v)  =  I  _ _ - 

'  7T  ,  ,2  2 

(v-v  )c  +  or 
'  o' 

and  the  full  Lorentz  line  shape  has  the  form 

Sa  Sa 


(22) 


(v-v0)2  +  a2 


(v+vq)2  +  a2 


(23) 


where  S  is  the  intensity,  a  is  the  halfwidth,  and  vQ  is  the  center  fre¬ 
quency  of  the  absorption  line.  For  infrared  frequencies,  v+vQ  >>  a,  and 
the  second  term  in  Equation  (23)  is  negligible. 

The  fact  that  the  experimentally  observed  absorption  is  greater 
than  the  calculated  absorption  is  just  the  opposite  from  the  continuum 
absorption  caused  by  CO,,  in  the  14  urn  region.  For  the  C02  continuum,  the 
absorption  coefficient  is  smaller  than  the  calculated  absorption  based 
on  wings  of  Lorentz  lines.  The  suggestion  has  been  made  that  the  increas¬ 
ed  absorption  for  the  case  of  water  vapor  is  due  to  absorption  by  (H^O)^ 
dimer  molecules  [17]. 

The  water  dimer  molecule  is  in  equilibrium  with  water  vapor  through 
the  dissociation-association  reactions  (H^O)^  +  2  H^O.  The  equilibrium 
constant  for  these  reactions  is 

r»2 


H20 


'Eq 


0 


(24) 


Comparison  of  experimental  values  of  C°(v)  at  296  K  with  calculated  valu 
continuum  absorption  is  due  to  the  extreme  wings  of  simple  Lorentz  lines 
entz  lines.  Reproduced  from  Reference  Ll5l. 
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where  the  subscript  D  refers  to  the  dimer. 

Energy  transmitted  by  a  gas  of  dimer  molecules  is  expected  to  fol¬ 
low  the  usual 


Tg(v)  =  exp(-kD(v)uD) 


(25) 


where  kp(v)  is  the  absorption  coefficient  of  the  dimer  gas  and  Uq  is  the 
optical  thickness  which  can  be  written  by  Equation  (5)  as  Up  =  Wq-L. 

From  Equation  (25),  kg(v)  can  be  written  as 


Mv)  =  TIT  ln  Tn(v)  • 


(26) 


If  the  water  dimer  is  indeed  causing  the  absorption  in  the  8  urn  to 
12  pm  region,  then  kg(\>)  is  the  experimentally  measured  absorption  coeffi¬ 
cient.  However,  this  coefficient  is  measured  experimentally  with  respect 
to  the  concentration  of  water  vapor  (uH  Q),  not  with  respect  to  dimer 
concentration.  From  Equation  (24),  Ug  can  be  converted  to  uH  g. 


ud  =  ir|-  uh2o  ph2o 


(27) 


Substitution  of  this  expression  into  Equation  (26)  gives 


k0(v)  =  tt  ln  Tn(v) 

h2o  h2o 

Rearranging  terms  gives 


-1  kn( v) 

>  r°<”)  “  77  V 


(28) 


(29) 


The  left  hand  side  of  this  equation  is  the  expression  which  would  be 
measured  as  in  Figure  1,  assuming  the  continuum  absorption  is  caused  by 
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water  vapor.  Ignoring  dimer-dimer  contributions  to  kg(v),  this  expres¬ 
sion  is  linear  in  the  water  vapor  pressure.  This  is  the  same  pressure 
dependence  as  Burch  used,  assuming  that  the  continuum  was  caused  by  ex¬ 
treme  wings  of  very  strong  lines.  Therefore,  a  study  of  the  pressure 
dependence  of  continuum  absorption  in  the  8  urn  to  12  urn  region  is  not 
expected  to  support  one  absorption  mechanism  over  the  other. 

The  temperature  dependence  of  the  continuum  absorption  has  been 
cited  by  many  as  strong  evidence  for  the  mechanism  of  dimer  absorption 
[17,  22].  If  dimers  are  responsible  for  the  absorption,  the  absorption 
coefficient  should  have  a  temperature  dependence  of  the  form  exp(E/RT), 
where  E  is  the  bonding  energy  of  the  dimer.  Results  of  temperature  stud¬ 
ies  by  several  workers  show  that  the  temperature  dependence  of  the  continu¬ 
um  absorption  coefficient  can  be  fitted  to  an  exponential  form.  Figure 
7  demonstrates  this  fit.  Table  3  shows  that  a  reasonable  estimate  for 
the  constant  E/R  is  approximately  1800  K.  On  the  basis  of  a  hydrogen- 
hydrogen  bond,  the  bonding  energy  for  the  water  dimer  has  been  estimated 
at  3  to  5  kcal/mole  Cl 7 ,  22].  Using  R=2  cal/mole°K,  these  bond  energies 
give  a  value  of  T  between  1500  K  and  2500  K. 

The  agreement  between  the  temperature  dependence  of  the  water  vapor 
continuum  and  the  temperature  dependence  of  dimer  formation  plus  the 
fact  that  continuum  absorption  is  much  stronger  than  predicted  from  Lor- 
entz  lines  can  be  taken  as  strong  evidence  in  favor  of  the  dimer  absorption 
mechanism.  However,  as  Burch  points  out  [26],  so  little  is  known  about 
the  true  shapes  of  the  extreme  wings  of  water  vapor  lines  that  their  con¬ 
tribution  cannot  be  ruled  out. 

Results  such  as  Equations  (9),  (17),  and  (19)  are  largely  empirical 
and  do  not  depend  on  a  particular  absorption  mechanism.  They  generate 
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curves  which  are  fitted  to  the  data  and  thus  represent  algorithms  for 
evaluating  the  absorption  coefficient  under  specific  environmental  condi¬ 
tions.  As  such,  these  results  do  not  possess  any  power  to  predict  re¬ 
sults  of  experiments  outside  the  limitations  of  the  variables.  Absolute 
determination  of  the  mechanism  of  continuum  absorption  in  the  8  urn  to 
12  urn  region  would  enable  predictions  of  the  behavior  of  this  absorption 
under  new  circumstances  to  be  made. 

One  idea  which  may  help  to  determine  the  absorption  mechanism  is 
to  study  samples  of  0^0  instead  of  1^0.  The  isotopic  shift  in  the 
band  from  6.3  urn  to  8.5  urn  would  effectively  bring  the  CO^  laser  emission 
frequencies  closer  to  the  band  center  for  a  different  look  at  absorption 
due  to  the  wings  of  strong  lines.  A  close  look  at  the  pressure  and  temp¬ 
erature  dependence  of  the  continuum  absorption  of  D^O  might  well  shed 
considerable  insight  on  the  continuum  absorption  of  H^O. 

The  Mater  Vapor  Continuum  in  the  Region  of  4  ym 

D.  Burch  [15]  has  studied  continuum  absorption  in  the  4  um  region 
using  a  long  path  absorption  cell  and  a  grating  spectrometer.  Continuum 
absorption  in  this  region  arises  from  two  effects.  One  is  water  vapor 
absorption  which  has  been  attributed  to  the  extreme  wings  of  strong  water 
lines  in  the  bands  at  6.3  urn  and  2.7  um.  The  other  is  pressure-induced 
nitrogen  absorption. 

Burch  found  the  experiment  difficult  to  perform  and  his  accuracy 
was  less  than  he  had  anticipated.  The  major  difficulty  was  due  to  the 
fact  that  the  absorption  in  this  region  due  to  the  water  vapor  continuum 
is  close  to  two  orders  of  magnitude  smaller  than  continuum  absorption  near 
10  \fn.  As  a  result,  he  had  to  use  very  high  temperatures  (up  to  428  K), 
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with  water  vapor  pressures  of  100  Torr.  This  led  to  problems  of  water 
vapor  absorption  onto  the  cooler  mirrors  which  reduced  reflectivity.  Fig¬ 
ure  9  is  a  plot  of  the  self-broadening  coefficient  C°(v)  for  water  vapor 
absorption  in  the  region  from  2400  cm’^  to  2800  cm'^  at  several  different 
temperatures.  The  magnitude  of  the  values  of  C°(v)  in  this  plot  should 
be  compared  to  the  magnitude  of  C°(v)  in  Figures  5  and  6. 

If  the  temperature  dependence  of  the  water  vapor  continuum  absorp¬ 
tion  in  the  4  urn  region  is  accurately  represented  by  Equation  (18),  then 
plots  of  ln(C°(v))  as  a  function  of  1/T  for  various  frequency  positions 
should  be  linear.  Figure  10  shows  such  a  plot  for  four  selected  frequen¬ 
cies.  The  linear  fits  to  the  data  are  very  good,  and  the  slopes  are 
approximately  equal,  indicating  that  the  TQ  value  in  this  region  is  not 
frequency  dependent.  The  lines  must  be  extrapolated  outward  to  find  C°(v) 
at  room  temperature,  since  this  quantity  cannot  be  measured  directly. 

Several  studies  [27,  28]  of  the  propagation  of  radiation  from  the 
deuterium  fluoride  laser  have  used  the  results  of  Burch.  F.S.  Mills  [27] 
has  stated  that,  on  the  basis  of  his  OF  laser  studies,  the  extrapolation 
done  by  Burch  was  probably  correct.  However,  Mills  recommended  that  this 
problem  should  be  studied  carefully  with  the  use  of  a  spectrophone. 
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Figure  10.  Semi- logar i thmic  plots  of  C°(v)  versus  1/T  for  four  differ. 

ent  frgquencies.  The  plotsSare  extrapolated  to  T=296  K  to 
find  Cj(v)  at  room  temperature. 
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SECTION  III 

THE  No-BROADENEO  WATER  VAPOR  ABSORPTION  LINE  SHAPE  AND 
INFRARED  CONTINUUM  ABSORPTION 

Part  I:  Theoretical  Development 

INTRODUCTION 

Strong  absorption  bands  of  water  vapor  and  carbon  dioxide  separate 
spectral  regions  of  high  infrared  transmi ttance  within  the  lower  at¬ 
mosphere  into  several  spectral  windows.  It  has  long  been  known  that 
attenuation  of  radiation  within  these  transmittance  windows  is  greater 
than  predicted  from  calculations  of  local  Lorentzian  absorption  line 
contributions  and  from  possible  aerosol  contr ibut ions.  The  residual 
absorption  in  excess  of  the  predicted  amount  has  been  called  continuum 
absorption  since  it  is  a  slowly  varying  function  of  frequency.  This 
continuum  absorption  represents  a  depression  in  the  spectral  background 
from  the  100%  transmittance  level,  and  has  been  linked  to  water  vapor 
content  along  the  absorbing  path. 

Continuum  absorption  in  the  8-14  um  region  has  been  the  topic 
of  a  great  deal  of  study  [1-5].  Investigations  of  its  behavior  as 
a  function  of  total  atmospheric  pressure,  partial  water  vapor  pressure, 
and  temperature  have  been  made.  Similar  measurements  have  been  per¬ 
formed  in  the  3.5-4  um  region  [6,7],  although  these  experiments  have 
been  more  difficult  because  the  continuum  absorption  in  this  region 
is  mucii  weaker  than  the  absorption  around  10  um. 

The  exact  mechanism  of  water  vapor  continuum  absorption  within 
a  transmittance  window  has  not  been  understood.  Originally  it  was 
thought  to  be  the  result  of  wings  of  very  strong  water  vapor  lines 
located  in  the  absorption  bands  bordering  the  transmittance  window. 


However,  it  was  shown  [2]  that  self-broadened  Lorentzian  absorption 
line  wings  were  not  adequate  to  model  the  10  N  continuum  absorption 
caused  by  pure  water  vapor  samples.  Furthermore,  absorption  in  the 
wings  of  Lorentzian  lines  near  10  um  increases  with  increasing  temperature 
while  experiments  show  the  measured  absorption  actually  decreases  with 
increasing  temperature  [8,93. 

Based  on  these  inadequacies  of  the  Lorentz  absorption  line  wing 
theory  to  model  the  water  vapor  continuum,  other  mechanisms  for  the 
absorption  were  sought.  Infrared  absorption  by  water  dimers  [ 1 0 , 1 1 3 
or  by  more  complex  molecular  clusters  [123  have  been  proposed  as  possible 
attenuation  mechanisms.  It  is  argued  that  water  vapor  aggregates  could 
have  unresolved  absorption  bands  throughout  the  infrared  spectral  region 
and  although  the  fractional  concentration  of  aggregates  would  be  small 
compared  to  the  concentration  of  H.,0  molecules,  the  absorption  could 
still  be  large  enough  to  be  detected  in  the  transmittance  windows. 

Perhaps  the  most  convincing  argument  in  favor  of  the  attenuation 
being  caused  by  water  vapor  aggregates  is  that  this  mechanism  predicts 
the  observed  negative  temperature  dependence.  It  is  argued  that  since 
the  attenuation  is  directly  proportional  to  the  number  of  aggregates 
present,  which  in  turn  is  proportional  to  exp(-AH/kT)  where  aH  is  the 
change  in  enthalpy  on  aggregate  formation  and  is  negative,  the  resulting 
attenuation  will  exhibit  a  negative  temperature  dependence.  Furthermore, 
measurements  of  the  temperature  dependence  of  water  vapor  continuum 
absorption  in  the  8-14  um  spectral  region  are  in  good  agreement  with 
predictions  based  on  hydrogen-hydrogen  bond  dissociation  rates  [133. 

Of  course  there  is  always  the  possibility  that  both  absorption 
mechanisms  are  occuring  simultaneously.  Measurements  of  the  temperature 


38 

dependence  of  continuum  absorption  at  1203  cm  *  by  Montgomery  [14  ] 
showed  that  as  the  temperature  of  the  water  vapor  sample  was  increased 
from  room  temperature,  the  absorption  reached  a  minimum.  This  behavior 
was  interpreted  by  Montgomery  as  showing  the  competitive  effects  of 
dimer  absorption  and  tar  wing  absorption.  He  suggested  that  initially 
the  measured  absorption  decreased  as  the  temperature  increased  because 
of  dimer  dissociation .  Above  400  K,  however,  with  dimers  depleted,  far 
wings  of  the  strong  absorption  lines  of  water  vapor  began  to  dominate 
and  push  the  absorption  up  slowly. 

Figure  1  demonstrates  that  the  accumulative  effects  caused  by 
the  overlap  of  absorption  line  wings  can  be  significant.  Three  spectra 
of  the  4.3  pm  absorption  region  of  the  COp  molecule  are  shown.  The 
background  spectrum  shows  the  absorption  caused  by  atmospheric  CO^ 
in  the  few  meters  of  path  between  the  infrared  source  and  the  entrance 
to  our  evacuated,  long-path  absorption  cell.  The  high-frequency  edge 
of  the  band  of  CO.,  shows  an  abrupt  cut-off  of  the  absorption  near 
2380  cm~*  which  is  characteristic  of  the  band  head. 

A  spectrum  recorded  with  80  torr  of  pure  C0^  in  a  300  m  path, 
however,  shows  that  the  absorption  caused  by  this  band  can  extend  beyond 
the  band  head.  Ignoring  the  very  weak  resonance  absorption  band  struc¬ 
tures  which  become  visible  at  this  optical  depth,  it  is  clear  that 
a  structureless  absorption,  manifested  as  a  slowly  varying  depression 
of  the  background  exists  even  past  2450  cm"*.  Filling  the  absorption 
cell  to  740  torr  with  pure  CO-,,  as  shown  in  the  third  spectrum,  causes 
the  influence  of  this  absorption  band  to  extend  past  2550  cm”*. 

This  increase  in  absorption  at  frequencies  above  the  band  head 
near  2380  cm"*  is  the  result  of  the  accumulation  of  far  wings  of  the 
absorption  lines  of  the  ^  band.  To  model  this  absorption  properly. 
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Figure  1.  Far  wing  absorption  by  CO ^  molecules. 


the  shape  of  the  absorption  line  profile  for  each  CC^  line  in  the  ab¬ 
sorption  band  is  needed.  Clearly,  an  incorrect  line  profile  which 
underemphasizes  each  contribution  to  the  total  absorption  in  the  far 
wings  would  result  in  excess  absorption  appearing  to  be  present  in 
the  measured  spectrum. 

The  spectrum  of  CO^  was  used  to  visualize  the  continuum-like 
absorption  which  can  be  caused  by  fay-  wings  of  resonant  absorption 
lines.  For  the  case  of  water  vapor,  it  is  difficult  to  find  clear 
spectral  regions  free  from  local  lines  in  which  to  observe  absorption 
line  wings.  Furthermore  the  physical  properties  of  water  vapor  prevent 
the  use  of  reasonably  large  pressures  (near  1  atm,  for  example)  without 
high  temperatures.  Thus,  it  is  difficult  at  tunes  to  separate  pressure 
dependent  and  temperature  dependent  effects  experimentally. 

In  order  to  model  water  vapor  continuum  absorption,  it  is  necessary 
to  determine  the  contributions  made  by  the  accumulation  of  the  far 
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wings  of  pressure-broaded  water  vapor  absorption  lines.  Only  then 

can  the  importance  of  absorption  mechanisms  based  on  molecular  aggregates 

be  established. 

In  this  study  we  have  begun  to  investigate  the  importance  of 
far  wing  absorption  for  characterizing  the  water  vapor  continuum. 

An  infrared  line  shape  for  pressure-broadened  water  absorption  has 
been  derived  which  is  valid  under  tropospheric  conditions.  It  was 
found  during  the  study  that  the  traditional  assumptions  which  lead 
to  the  Lorentz  line  shape  could  not  be  used  when  considering  the  far 
wings  of  water  vapor  absorption  lines. 

The  total,  normalized  line  profile  was  used  to  predict  the  ab¬ 
sorption  at  specific  frequencies  within  the  8-14  urn  and  3.8-4  urn  trans¬ 
mittance  windows.  The  experimental  methods  and  the  results  are  dis¬ 
cussed  in  detail  in  anolf.i  r  report  Ci53.  -in  this  section  we  shall  discus: 
the  development  of  the  infrared  absorption  line  profile  for  water  vapor 
molecules. 

THEORY 

The  absorption  coefficient,  k(v),  for  attenuation  caused  by  ab¬ 
sorption  lines  can  be  written 

k(  v)  =  l  si J'i(  v  v0i )  d) 

i 

where  is  the  line  strength  of  the  i ^  absorption  line,  and  j.(v- 
v  ■)  is  the  line  shape  function  of  the  ilh  absorption  line  centered 
at  frequency  v  It  can  be  assumed  that  any  differences  in  the  line 
shapes  of  different  water  vapor  lines  can  be  expressed  parametrically 
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within  the  same  overall  function.  That  is,  J-j(v_v0i)  Is  better  written 
as  j ( v  —v o ^ ^ -j )  where  P-  represents  a  set  of  parameters  which  represent 
the  ith  absorption  line. 

The  line  shape  function  is  normalized  by  the  requirement 

oo 

j ( v- x'oi  ;p  ■  )dv  =  1  (2) 

Jo 

for  all  i. 

Any  description  of  the  absorption  coefficient,  k(v),  then  can 
be  divided  into  a  line  strength  term  and  a  line  shape  term. 

The  Hamiltonian 

The  starting  point  for  the  development  of  the  spectral  line  profile 
is  with  the  total  Hamiltonian  describing  the  molecular  system.  In 
this  case,  the  system  consists  of  a  number  ng  of  gas-phase,  interacting 
molecules  which  can  absorb  infrared  radiation.  Also  present  are  other 
molecules  which  can  interact  with  the  absorbing  molecules  via  collisions. 
Thus,  the  Hamiltonian  can  be  written  as 

HTot(0  *  %  ‘Vp*  Vl)  +  H‘(t)  (3) 

where  H,  is  the  Hamiltonian  for  the  absorbing  molecule 
dl 

H.  is  the  Hamiltonian  for  the  buffer  molecules 
b 

Hp  is  the  Hamiltonian  for  the  photon  field 

H  (t)  is  the  interaction  Hamiltonian  between  the  absorbing 
djP 

molecule  and  photon  field 
Hc(t)  is  the  collision  Hamiltonian  which  describes 
molecule-molecule  interactions. 
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Time  dependent  terms  are  explicitly  indicated.  The  collision  Hamiltonian 
will  be  written  to  account  for  absorber-absorber  molecule  interactions 
(spectral  self-broadening)  and  to  account  for  absorber-buffer  molecule 
interactions  (spectral  foreign-broadening). 


indicates  the  specific  molecule  considered  to  be  doing  the  absorption 
while  the  n  -1  identical  molecules  are  considered  for  the  moment  to 

a 

be  perturbers.  The  unperturbed  Hamiltonian  for  the  absorbing  molecule 
is 


=  H. 


+  H. 


+  H. 


+  H. 


ROT 


VI B 


EL 


TRN 


where 


(4) 


H  describes  the  rotation  energy  of  the  molecule 
aR0T 

H  describes  the  vibrational  energy  of  the  molecule 
vIB 

H  describes  the  electronic  energy  of  the  molecule 
dEL 

and  H  describes  the  translational  energy  of  the  molecule. 

“TRN 

The  first  three  terms  are  fundamental  to  the  study  of  molecular  spec¬ 
troscopy  and  to  the  production  of  spectral  lines.  The  term  H  will 

aTRN 

not  be  directly  included  in  this  analysis.  Motion  perturbed  line  pro¬ 
files  such  as  the  Doppler  profile  can  be  derived  from  the  term. 

In  the  total  Hamiltonian,  the  term  is  a  many-particle  Hamiltonian. 
This  includes  foreign  molecules  and  molecules  identical  to  the  absorber. 

It  can  be  written 


Hb  =  l  Ha  + 

b  j=2  aj 


<  n^ 

11%. 

k=l  m=l  Dkm 


(5) 


Here  n.  is  the  number  of  absorbing  molecules,  k  represents  the  number 

u 
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of  types  of  foreign  molecules,  and  n^  is  the  number  of  molecules  of 
the  kth  type.  Each  Hamiltonian  is  structured  in  the  same  way  as  Equation 
(4);  that  is  the  rotational,  vibrational,  electronic,  and  translational 
energies  are  represented. 

Hp  is  the  free  space  photon  Hamiltonian,  written  as 


H  =  l  Mq)  [a*(q)a  (q)  +  l2  ]  (6) 

H  0  0 

qo 

where  ( q )  is  the  photon  creation  operator  and  a,,(q)  is  the  photon 
annihilation  operator  for  a  photon  of  polarization  ’  and  wave  number 
q  corresponding  to  the  angular  frequency  w=  ti|q|. 


The  interaction  Hamiltonian,  H  ft), 

a^p 

and  the  photon  field  can  be  written  [16] 


between  the  absorbing  molecule 


Ha  (t)  =  ^-[  P  •  A(r,t)  +  A(r,t)  •  p] 

0^  J.  _v 

+  Z— ?  A(r,t)  .  A(r,t)  (7) 

2mc^ 

wherelp  is  the  particle  momentum  operator,  and  A("r, t)  is  the  photon 
field  magnetic  vector  potential  operator,  written  as 


A(r,t)  = 


2jHic  1 2 
qV 


a0(q)>,(q)exp(i{q  r-  ft)) 


+  an  (q)>  *(q)exp(-i  (q  r-'-'t) ) 


-} 


(8) 


In  this  expression  T(q)  is  a  unit  vector  perpendicular  to  the  free 


space  photon  propagation  direction  and  represents  the  polarization 
of  the  electromagnetic  radiation. 


To  first  order,  the  interaction  Hamiltonian  becomes 


H,  (t)  =  Hab  exp(-iwt)  +  H^m  exp(iwt) 
a^p  a^p  a^p 


(9) 


where  the  A  •  A  term  has  been  dropped.  The  term  Hab  represents  the 

a^p 

collection  of  operators  and  factors  which  multiply  the  negative  ex¬ 
ponential  after  Equations  (7)  and  (8)  have  been  combined.  This  term 

represents  the  absorption  of  radiation  by  the  molecule  [17].  Similarly, 
pm 

the  term  H  represents  the  emission  of  radiation  by  the  molecule, 
a^p 

Finally,  Hc(t)  represents  the  intermolecul ar  potential  which 
couples  the  absorbing  molecule  to  molecules  surrounding  it.  It  is  the 
least  understood  term  in  the  total  Hamiltonian,  and  yet,  it  is  the 
most  important  when  attempting  to  understand  the  nature  of  resonant 
line  profiles.  A  general  many-body  expression  can  be  written  as 


Hc(t)  ■  1 


j=2  Vl 


k  n. 

(ti  *  i  l 


H.  (t) 
k=l  m=l  Dmkal 


(10) 


when  restricting  the  discussion  to  binary  collisions  only.  Again, 
k  is  the  number  of  different  foreign  gases,  and  n^  is  the  number  of 
molecules  of  the  kth  gas.  Figure  2  illustrates  the  geometry  involved 
in  such  an  interaction.  For  electrically  neutral  molecules  (no  ionic 
interaction)  with  permanent  electrostatic  potentials,  a  single  Hamiltonian 
in  Equation  (10)  takes  the  form  [18] 
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rmk,rl^  ^Electrostatic  +  ^Magnetostatic  >  Rot 

^Dispersive  ’^Inductive  +^repulsive 

As  shown  in  Figure  ?,  R(t)  is  the  interrno locular  distance  between  the 
absorber  and  perturbor.  The  quantities  and  r^  locate  the  charged 
particles  within  the  molecular  coordinates  of  the  foreign  and  absorber 
molecules  respectively.  The  magnetostatic  potential  function  is  sig¬ 
nificantly  smaller  than  the  elect  rost.at  ic  potential  and  this  i*  ignored 
hereafter. 
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A  brief  description  of  the  potential  terms  is  necessary.  A  multi¬ 
pole  expansion  of  the  electrostatic  potential  is 


Electrostatic 


yVl j  +  WFmk»rl>  +  WFmk’rl>  +  WV.x’V 


R3(t) 


R4(  t ) 


R5(t) 

(12) 


where 

^DD^mk’^l^  is  the  dipole-dipole  interaction  term 
VDq(rrn|<,ri)  is  the  dipole-quadrupole  interaction  term 
^Do(r?mk,r'l)  is  the  dipole-octopole  interaction  term 
and  vqq( ^mk ’ ^1 )  is  t,ie  Quadrupole-quadrupole  interaction  term. 

These  potentials  can  be  expressed  as 

VD[fmk’Fl)  =  umkl  lf(  mk’Vmk’7) 

VDQ(rmk’rl)  =  l,mk^l9l(  mk’";l’rmk’r^  +  ^mkl'l92(Rmk’ni>rmk’r) 
VDD(rmk’rl)  =  llmk^lhlFmk’nrrmkr)  +  ^mklllh2(f!mk’fil’rmk’r) 
and  =  Wlq("mk-l’7mk’7)- 

The  is  the  dipole  moment  of  the  mth  molecule  of  the  kth  type,  and  is 

the  dipole  moment  of  the  absorber.  Employing  the  same  convention  with 
the  subscripts,  $  and  IT  represent  the  mean  quadrupole  and  octopole 
moments  of  the  molecules,  respectively.  The  terms  and  ^  represent 
the  parameters  necessary  to  establish  the  orientation  of  the  multipole 
moment  of  both  the  perturber  and  absorber  molecules  in  a  space-fixed 
coordinate  system.  These  parameters  are  arguments  of  expressions  f, 
g,  h,  and  q  which  determine  the  angular  dependence  of  the  electrostatic 
interaction  potentials. 


47 


The  molecules  are  rotating  during  the  collision.  The  total  ro¬ 
tation  during  interaction  will  depend  on  the  rotational  energy  of  the 
molecule  and  the  duration  of  the  collision.  The  forces  produced  by 
the  electrostatic  potential  will  attempt  to  orient  the  molecules  in 
a  minimum  potential  energy  conf iguration.  A  weighted  average  over 
the  rotation  is  performed  in  order  to  account  for  the  ability  of  the 
absorber  molecule  to  achieve  this  preferred  orientation.  The  average 
is  denoted  by  <  >^Qt,  and  is  written 


<V 


E^Rot  = 


Jd,  VE,e»p(-VE ,/H)p(--0.td.4.f” ) 
fd>:  0xp(-VE]/kT)p(i:o,td,f£k,flJ) 


(13) 


Here  dil  is  the  differential  element  of  solid  angle  for  all  possible 

orientations  and  p (  — 0 , , f ^ » f i )  is  a  pulse  function  with  unit  amplitude 

and  a  pulse  duration  whicn  depends  explictly  on  the  initial  orientation 

n  ,  the  collision  duration  time  td,  and  the  rotational  frequencies 

R  R 

of  the  perturbing  and  absorbing  molecules  f(J)k  and  f^,  respectively. 

For  collisions  with  large  impact  parameters  the  collision  duration 
time  is  very  small,  and 

‘'TiWv  on 


For  collisions  with  small  impact  parameters  the  electrostatic  potentials 
can  be  so  strong  that  the  interaction  can  no  longer  be  considered  to 
be  a  small  perturbation  to  the  unperturbed  Hamiltonian  of  the  absorber 
molecule.  This  occurs  during  close  encounters  and  during  the1  formation 
of  liquids  and  solids.  The  ability  of  the  molecule  to  rotate  is  in¬ 
hibited  and  the  absorber  molecule  takes  a  relatively  fixed  orientation. 
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At  this  extreme,  the  average  over  orientations  breaks  down  and  in  fact 
the  formulation  of  the  unperturbed  Hamiltanian  as  Equation  (4)  is  no 
longer  valid. 

Other  long  range  electrostatic  potentials  can  be  of  interest 
to  molecular  collision  dynamics.  The  most  important  of  these  are  in- 
duction  and  dispersion  potentials,  both  of  which  vary  as  R~°(t).  In¬ 
duction  potentials  represent  the  interaction  of  a  permanent  electric 
multipole  of  one  molecule  with  an  induced  multipole  of  another  molecule. 
The  resulting  potential  is  attractive.  Dispersion  potentials  arise 
from  the  random  motion,  within  the  bounds  of  the  Uncertainty  Principle, 
of  electrons  about  the  nuclei  generating  instantaneous  multipoles. 

A  discussion  of  these  potentials  is  given  elsewhere  Cl8  ]  and  approximate 
forms  are  [19  ] 


V 


Inductive 


u2,«.  1 
mk  1 


R6(t) 


(15a) 


V 


Dispersive 


-22 
ae  <r  > 

2R6( t) 


(15b) 


where  a  is  the  mean  electric  polarizabi 1 i ty  of  the  molecule,  and  e 
is  the  charge  on  an  electron. 

The  final  term  in  the  collision  Hamiltonian  as  written  in  Equation 
(11)  is  Vpepu]sjVP-  During  very  close  encounters  between  molecules, 
repulsive  forces  such  as  electron  cloud  interaction  can  occur.  For 
this  study,  such  potentials  were  not  explicitly  considered. 
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THE  TRANSITION  RATE 

Following  the  work  of  Anderson  [  20 J  and  of  Tsao  arid  Curnutte 
[213,  the  rate  for  a  radiation  absorption  transition  from  a  lower 
molecular  energy  level  to  an  upper  molecular  energy  level  separatee 
by  Eu-E  1  =  is 


t-u  =  ^  S‘-1  J'0dI  e»p(-itoT)Ir[ooH^pS;1(  t)l^pSc(e]l,  (16) 


where  S  ( t )  is  the  time  development  operator  for  the  collision  Hamiltonian 
at  time  T  and  is  written 


Sc( t)  =  exp  ~  /  Hc(x)dT. 

0 

In  Equation  (16),  p0  is  tlie  density  matrix  generated  by  the  initial 
states,  and  Aui  =  wai0  is  the  amount  of  off-resonance  between  tfie  radi¬ 
ation  frequency  ;o  and  energy  level  spacing,  This  equation  can  be  written 
in  a  compact  form  as 


r 

H-hj 


— ~  Re  {  J  c|x  exp(-iA.i>T)  * 

W  0 


r,v.vi 

U 


1  ,  ,  l'  |  | <  1  > /  ,  M  w  . 

<d«,1  Vv^'qult  '  ijp  du',u'usu  q, 


‘VuVu'V.i 
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In  this  expression, 

g^  is  the  degeneracy  of  the  lower  level, 


g  -  pVibpRot 


is  the  density  matrix  over  vibrational  and 
rotational  states 


pp  is  the  density  matrix  of  the  photon  field 

nqair  nqat  are  Ptnoton  numbers  when  the  system  is  in  the  upper 

(final)  and  lower  (initial)  states  respectively, 

and  p  ^  is  the  density  matrix  over  vibrational  and  rotational  states 
s  n 

of  the  broadening  molecule. 

The  “a1'  quantum  numbers  represent  all  quantum  numbers  other  than  J 
and  M  necessary  to  describe  the  upper  and  lower  levels  of  the  absorbing 
molecule.  The  terms  s^  and  represent  all  quantum  numbers  necessary 
to  describe  the  lower  and  upper  levels  of  the  broadening  molecule. 

The  summations  in  Equation  (17)  are  over  states  and  should  be  written 
as 

l  -  l 
1  V*Mi 


v°t 

where  q  ^  is  the  initial  wavenumber  state  of  the  photon  and  is  the 
polarization  state  of  the  photon.  A  similar  expansion  can  be  written 
for  the  upper,  final  state. 

Abreviating  the  notation,  let  ?.  =  a  J  ^  ^  and  u  au'*u^u'  Realizing 
that  H^p  is  independent  of  and  sj^.  Equation  (17)  can  be  written 
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?  faVJ.M.  f  I  J  W  ?  .  *4  4*t*s 


g  a  L  p . qo* '  alP'u  "qau 

\iV'  I  q.o.  1' 


a‘ jV  i  V°u 
u  u  u 


X  j  d,  e-’A“’<u  SJ  nqou|S-‘(,)H^p  Sc(T)'t^  nq3l>  (18) 


Reducing  the  sums  for  nontrivial  matrix  elements 


r  £  71  r> 

« =  7 Re 


[  p,|-  I  l  p  <?n  ]Hf  jun  >  y  rs‘ 

[a  J  1  9  it  M  M  q  a,  p  qolV  qcu  Lh 


V  l  '  V  u  ^ 

q  ,  a 
Mu*  u 


-  f  dx  e'lAwT  < u  sbn  ! S”1 ( T ) Hab  Sr(T)|tsbn  > 
tt  J  Q  qcu  c  v  '  a^p  c” ' 1  qaf. 


The  subscripts  u  and  i  have  been  dropped  from  the  SD  since  the  Kronecker- 

5  in  Equation  (18)  requires  J-=u.  By  inserting  unit  operators  between 

the  time  development  operators  and  the  interaction  Hamiltonian  Hab 

a^p 

the  second  matrix  element  expression  in  Equation  (19)  becomes 


^"qaulSc’ta^p  Sc(t)  |e\ot>- 

I  <  A,uM')|u'Sb'n'qou’-J'sb' 


«b\-. 


t,sh"  "qrt,Sc< 
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Since  S  ( t )  is  independent  of  the  photon  field,  and  Hdb  does  not 
operate  on  states  of  the  broadening  molecule,  then  sb'  =  sb",  n 


qo2, 

nqjji’  anq  nqou  ~  nqou* 

The  adiabatic  assumption  [21,22]  is  now  made  requiring  u  =  u‘, 
b  b 1 

s  -  s  ,  and  1=1'.  By  this  approximation,  the  right  hand  side  of 


Equation  (20)  becomes 

<USb|  S‘l(T)|  UsbxUnqou|  Hbbp|  inqolxlsb|  SC(T)|  lsb>. 


(21) 


Substituting  this  result  into  Equation  (19),  the  transition  rate  can 
be  written 


r  =  2^ 
■  *2 


£ 


i 

9 


4-  l 


1M»Mu 

V°u 


l 


p  |  <un  j  Hab  I  Jtn  _>|^  — 
p1  qou1  a,  p1  qoi.  1  it 


'od 

dr 

o 


X  Psb  <-b|s-1(T)|usb><Ub|Sc(T)Usb>.  (22) 

Assuming  binary  collisions  only,  the  matrix  elements  of  the  time  develop¬ 
ment  operator  Sc(x)  can  be  written 


-bi  c  2_\in.b.  _  ^fflr.aic  i  „  _a.  na  bl, 


bl  nl 


Cits  |  Sc(x)|  fts  >  =  <ls  I  scaial  g-s3>  <£s  lsr<,  h.l£s  A>  ••• 


caabl 


_ <£s  k| S-  .  |e  bk>nk 

1  calbk'  s 


(23a) 


~bi  c-1 


<usu|  S‘x(t)|  usu>  =  < u s a j  S ” ^  lusa>n<1  '  <us  ^S'1  .  us  l> 


bl  nl 


C  3  j  J 


CJabl 


bki  -1  bk 

-..<9.1.  S .[  ,  )  Sts  k> 

'  ca, n  ' 

1  k 


-  f  Aojt 


(23b) 
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In  these  expressions,  the  terms  similar  to  S  ,  represent  quantities 

c“lbk 

averaged  over  impact  parameters  and  summed  over  collisional  events. 

The  sums  over  broadening  states  sb  in  Equation  {22)  can  be  expanded 


l  p  b<usb|S^T)|usb><i,sb|  Sc(T)jHsb  >- 

sb  S 

-  n. 

I  Ds3«usaj  .(t)  |usa><lsa|  Sca  a(T)l  fl'Sa>>R  ,  v.t  ‘ 

a  A  X  L)  U  _ 

Is 


xr  I  0  bk  1  bk  bk  ,  bk  1  k 
1  AsV*  'Vv.tl 


The  symbol  <  x.  indicates  the  average  over  impact  parameters  and 

Ro’v’to 

velocities  and  a  sum  over  collision  events.  It  is  written  mathematically 
for  arbitrary  A  as 

co 

<»* ■  v  «*1/2  ( *  >3/2  f  R  «o4l 

o’  ’  o  Jo  M in  b 


where  V  is  volume,  u  is  the  reduced  mass  of  the  binary  system,  T  is 
temperature  and  k  is  the  Boltzmann  constant.  Also,  because  na>>l  the 
approximation  na-l~na  was  used  in  Equation  (24). 

It  is  convenient  at  this  point  to  define 

la  1  I 


Ca1bk(T)  =L  ^  ScLlJT*,s  k><5-S  k!Sca.h.  k>^0>v,to. 


Ik  \  s 
s 


CW 


C(i,h, 
1  k 
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Then 


Re  I  dr  e'iAa)TC(i)  =  <5(Aw), 

71  Jo 

and  the  transition  rate  as  expressed  in  Equation  (28)  takes  the  familiar 
form  of  Fermi's  Golden  Rule. 

The  right  hand  side  of  Equation  (28)  is  the  product  of  two  terms. 

The  first  term  is  interpreted  as  the  transition  strength,  while  the 
second  term  is  an  off-resonance  amplitude  correction  to  the  transition 
strength.  This  second  term  then  is  the  frequency  profile  of  the  transition, 
otherwise  referred  to  as  the  lino  shape.  We  write 
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j|?(  Ato)  =  -  Re 

C  7T 


d  t  exp(-i  Ao)t)C(  t)  , 


(29) 


where  the  subscript  c  refers  to  the  collision  broadened  line  shape. 


THE  CORRELATION  FUNCTION 

C(  t)  as  defined  in  Equation  (27)  is  interpreted  as  the  autocor¬ 
relation  function  between  the  state  of  the  molecular  system  at  time 
t=0  and  time  t=T.  Just  as  an  interferogram  produced  by  a  Michelson 
interferometer  is  the  autocorrelation  function  of  the  incident  electric 
field  and  is  symmetric  about  the  zero  retardation  point,  so  the  function 
C(x)  is  symmetric  about  r=0.  Using  this  symmetry,  Equation  (29)  can 
be  written 

R  Dp  f® 

Jc(Aw)  =  Ytt  dT  exp(-iojr)exp(im0T)C(  t) 

'  —  cn 


=  Re  f^'1  (exp(  i uiot)C(  t)  } 


(30) 


That  is,  the  line  shape  function  j^(Aw)  and  the  function  exp(iu  t) 

C(t)  are  related  by  the  Fourier  transform.  To  generalize  the  concept 
further,  let  jc(Aa>)  be  the  Hilbert  transform  companion  to  j^(Am). 

We  define 

j  (Au>)  =  «31(CxP(i^T(C(x))  =  j^(M  +  i  j*(Aw).  (31) 

and  j^(Aui)  represents  refractive  index  effects  caused  by  the  collision 
processes.  Since  the  goal  of  this  study  is  a  more  thorough  understanding 
of  the  water  vapor  absorption  line  shape,  the  retractive'  index  terms 
will  be  ignored  and  the  superscript  R  suppressed.  Thus 


jc(M  =  Re{f'1(exp(  iw0t)C(  i)) }. 


(3?) 
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This  equation  indicates  that  a  study  of  the  absorption  line  shape  depends 
on  an  understanding  of  the  correlation  function  C(t).  An  intuitive 
understanding  of  the  relation  between  the  spectral  line  shape  and  the 
correlation  function  is  shown  in  Figure  3.  Long  correlation  times 
t  imply  that  the  binary  collision  perturbation  was  weak  and  did  not 
greatly  upset  the  energy  levels  of  the  absorber.  Thus,  appreciable 
absorption  caused  by  a  transition  from  the  lower  energy  level  to  the 
upper  energy  level  will  occur  at  frequencies  co  close  to  the  unperturbed 
center  frequency  wQ.  In  the  spectrum,  the  absorption  line  is  narrow. 

If,  however,  the  correlation  time  for  the  energy  levels  is  short, 
a  strong  interaction  between  absorber  and  perturber  molecule  occurred. 
This  can  cause  the  absorber  molecule  to  absorb  radiation  at  frequencies 
to  which  are  far  from  wQ.  The  spectral  line  appears  broad. 

Of  course,  in  an  experiment  all  ranges  of  correlation  times  are 
experienced  by  the  ensemble  of  absorber  molecules.  The  total  line 
shape  is  a  composite  resulting  from  the  different  ranges  of  correlation 
functions. 

The  Lorentz  Line  Shape 

It  is  reasonable,  but  not  entirely  correct,  to  assume  that  the 
correlation  function  is  a  monotonical ly  decreasing  function  of  positive 
t.  That  is,  the  molecular  system  tends  to  become  less  correlated  as 
time  proceeds.  If  this  monotonic  decay  can  be  written  as  a  simple 
exponential  of  the  form 

C(t )  =  exp( -tr )  (33) 

Then  from  Equation  (3^),  the  line  shape  becomes 


jc(^)  = 


1 

IT 


d 

Auf’  +  o? 


(34) 
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Figure  3.  Comparison  among  the  intermolecular  potential,  the  correla¬ 
tion  function,  and  the  spectral  line  shape. 

which  is  the  familiar  Lorentz  line  shape  with  a  halfwidth  a.  The  Lorentz 
line  shape  is  an  adequate  model  for  the  absorption  line  profile  only 
when  the  correlation  function  can  be  expressed  as  an  exponentially 
decaying  function.  It  will  be  shown  later  that  this  assumption  can 
be  applied  to  the  water  vapor  line  shape  only  under  certain  conditions. 

For  the  most  part,  the  profile  of  water  vapor  absorption  lines  is  cecicedly 
non-Lorentzian. 

Col  1 ision  Dynamics 

For  this  study,  it  was  decided  that  the  best  approach  for  establishing 
an  absorption  line  profile  was  to  classify  the  correlation  functions 
into  three  general  categories  corresponding  to  short,  intermediate, 
and  long  correlation  times.  By  a  previous  discussion,  this  classifi¬ 
cation  divides  the  absorption  line  profile  into  three  regimes  corresponding 
to  far  wing,  intermediate,  and  near  line  center  domains.  Figure  4 
illustrates  the  division  of  the  collision  dynamics  implied  by  this 
approach.  Trajectories  labeled  1  and  ?  represent  possible  colltsional 
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© 


Figure  4.  Possible  collision  trajectories. 


paths.  Table  I  lists  the  general  characteristics  of  collisions  in 

each  region.  These  conditions  must  be  applied  to  the  correlation  function 
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TABLE  1 

CHARACTERISTICS  OF  COLLISION  DYNAMICS 


:  t  t  rj(7  ~v  • ;  I -• ;  .  • 

*v|.  /  r  /  r  r.’/'f  ;  <i/7  :>/ ■  >. ?; i  n<;  ■.  /  r  j  n-  •<  ;  :j.  .*,r <.-•< 

•  •  Biot  i'*n  o  t  <-<>}  1  i  r.  i,  >n 

i  *  .'•«*{  u*»  » -n  i  *  ‘  l  1  l i'lnr.  :  ■  rif-or-l  to 

U  /.it  ■  iJl 1  i  >  /  <  ■<  i J  ;  S  i  <  »ns 

i.  fc.’i.n  -;  .  t>!  I  '.r.umr, 

■l .  ■/.*>/  or !  •  !  it  :  t  )i,  t  i  pi,  > 

/'  f**t.  7/Jf  • -f  .*Ht  •*//.<  f « •  r»  •*;;'«  */j' j  u.'r  r  . 

.  Msiriu  i'oJcj  :  i  >}-h  in 

P  • ' *  <  •*  *  •'*’ '  *  •  •  <  *  -i ■  •’  :  /;■ !  j i  w  /;•. j 

J  .  BJlltij  i;  *.<  > ;  s  j ,  : j 

.  Pur. 4  f  i  o:j  '  *  •  >;j  '  ;L 

.*'hor f  . \ ■<:  i*  Lit  ion  ti no 


jc(Aw)  =  Re^'^expdajQtjCpyCT))  +  2_1(exp(  1u>0t)Ciw(t)) 
+  ^_1(exp(ioj0T)CNLC(T))} 


jcFW^)  +  +  ^cNLC^^  * 


(35) 


The  term  j  fl,(Auj)  describes  the  far  wing  collision  line  shape  and  is 


derived  from  correlation  functions  with  short  correlation  times.  The 


contributions  to  be  near-line-center  shape,  Jcnlc(^)>  come  ^rcnTI 
collisional  processes,  but  are  dominated  by  correlation  functions  with 
long  correlation  times.  The  intermediate  wing  region  Jcj^(Am)  represents 
a  transition  region  where  correlation  functions  with  both  long  and 
short  correlation  times  compete  equally.  Rather  than  solving  for  the 
line  shape  function  in  this  region,  boundary  conditions  that  the  total 
line  shape  function  and  the  derivative  of  the  function  must  be  continuous 
across  the  domains  will  be  applied. 
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Solutions  to  the  near- 1 ine-center  problem  and  the  far  wing  problem 
can  be  formulated.  The  intermediate-wing  region,  however,  is  not  solvable 
by  direct  methods. 

The  collision  Hamiltonians  for  H^O^O  and  H^O-broadener  binary 
interactions  become 


m  _  ,  ^iW^i  qi92  ,  , 

aal  R3(t)  R4  ( t )  Rot 

<rp 

H  a  +  _  2  4  "2  +  He2<r2> 

bai  <  R4( t)  ^ot  R6(t) 


(36a) 


(36b) 


The  terms  have  been  described  earlier  and  the  subscripts  1  and  2  repre¬ 
sent  water  vapor  and  broadener  molecule,  respectively. 

The  short  collision  times  characteristic  of  interruption  broadening 
in  Table  I  allow  the  approximation  of  Equation  (14)  to  be  made  for 
the  near  line  center  problem.  Retaining  only  the  leading  terms  in 
each  Hamiltonian  above  is  equivalent  to  assuming  that  the  H?0-H„0  inter¬ 
action  is  a  dipole-dipole  interaction,  while  the  H^O-N^  interaction  is  a 
dipole-quadrupole  interaction.  The  limiting  case  cf  weak  collisions  can 
be  modeled  with  the  Hamiltonians. 


aa. 


ba. 


R3  ( t ) 

(37a) 

uiQ2 

R4(t)  q 

(37b) 
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Strong  collisions  which  contribute  to  the  far  wings  of  absorption 
lines  may  or  may  not  require  evaluation  of  the  rotational  average  in 
Equations  (36a)  and  (36b).  The  importance  of  the  rotational  average 
can  be  evaluated  by  examining  the  collision  time  and  the  rotational 
time  for  the  molecules.  Choosing  the  broadening  molecules  to  be 
an  estimate  for  the  collision  duration  can  be  made,  based  on  the  average 
thermal  velocities  and  impact  parameters  of  strong  collision'  for  both 
H^O-HpO  and  H„0-Np  binary  collisions.  The  results  show  the  collision 
duration  to  be 


Using  a  rigid  rotor  picture  of  the  molecules,  typical  rotation  periods 
for  the  Np  molecule  are  between  10'^  and  10”^  sec,  while  for  the 
H^O  molecule  the  rotation  periods  are  between  10  and  10  sec. 

These  values  can  be  derived  from  a  purely  classical  view  of  the  molecule 
using  angular  momentum  values  which  are  derived  from  representati ve 
values  of  the  J  quantum  number  (J  £  5  to  20). 

Since  the  collision  times  for  HpO-Np  are  for  the  most  part  short 
compared  to  the  rotational  period  of  the  N?  molecule,  a  "snapshot"  view 
of  the  interaction  can  be  adopted.  That  is  the  Np  broadener  molecule 
does  not  rotate  significantly  during  the  interaction  and  therefore 
the  rotational  average  is  no'  performed.  Similar  results  are  derived 
when  considering  11^0-0^  collisions.  However,  for  HpO-H^O  collisions. 
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the  colliding  molecule  can  rotate  many  times  during  the  collision, 
makinq  the  rotational  average  important.  For  far  wing  absorption  line 
shapes,  the  binary  interaction  Hamiltonians  become 


H 


aa, 


2 

V(t) 


f> 

Rot 


2„fo  *  'i1e?<r2> 

R6(  t ) 


(38a) 


H  =  M1Q2 
bal  R4( t) 


(38b) 


The  average  over  rotation  as  defined  in  Equation  (13)  can  be  performed 
by  recognizing  the  analogy  between  this  and  orientation  polarizabi lity 
calculations  [23] . 


2 

hi 


2 

HI 


R3  ( t ) 


r>Rot 


R3(  t) 


L(x) 


(39) 


where  L(x)  is  the  Langevin  function  and  x  =  2p3/UTR3(  t)  ] ,  where  T  is  temper¬ 
ature.  F.gure  5  shows  a  plot  of  this  function.  If  x  is  small,  L(x)  ^  x/3, 
and  the  rotational  average  becomes 


2 

X 

R3( t)  >Rot 


3kT 


R6(  t) 


(40) 


which  is  similar  to  results  given  by  Eisenberg  and  Kauzmann  Cl9  ]- 

3 

The  next  term  in  the  expansion  of  L(x)  is  of  order  x  .  This  produces 

_  i  o 

a  term  in  the  Hamiltonian  of  order  R"  (t).  Terms  of  this  order  become 
important  during  very  close  collisions,  but  their  influence  drops  rapidly 
as  R(t)  increases.  For  this  study,  such  terms  are  not  being  considered. 
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Figure  5.  Langevin  function. 


Thus,  from  Equation  (38a),  the  self-broadening  influence  in  the 
far  wings  of  water  vapor  absorption  lines  is  controlled  to  first  order 
by  a  R  ^(t)  potential.  This  differs  from  the  self-broadened  nature 
of  the  near-1 ine-center  Hamiltonian  which  varied  as  R‘3(t)  for  large 
R(t).  The  forei gn- broadening  potent i al s  for  near- 1 ine-center  and  far 
wing  cases  both  vary  as  R’^(t). 
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(43) 
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To  obtain  separate  averages,  it  was  assumed  that  events  at  time  t 

and  time  x+dx  are  independent  [?1"|.  Within  the  limits  of  the  impact 

approximation,  the  second  average  is  independent  of  the  sum  over  quantum 

numbers  of  the  broadening  molecule  and  in  fact  is  equal  to  the  correlation 

function  C  ,  (t  )  at  time  x-  Thus,  Equation  (43)  can  be  written 
a^D 


C'aib^T+dx^  =  ^a^b^  Mx+x+dx) 


where  the  function  <f)(x-»x+dx)  simplifies  the  notation  of  the  remaining 
terms  in  Equation  (43),  and  is  dependent  only  on  dx.  A  solution  to 
Equation  (44)  can  be  written 


caib(,)  ”  Ca1b<0)<!’<P(-''«1bT) 


where  y,  .  can  be  shown  to  be 

D 


Y,  h  =  <1  -  y  p  b<ush|S’^  ,(  x-»x+dx)  just5><£sb|S  ,(  x+x+dx)  |  Hsn> 


b  „„b 


The  initial  condition  C  ,{0)  =  1  can  be  applied,  and  the  total  cor- 

a^b 

relation  function  for  a  binary  mixture  can  be  written  by  Equation  {??) 


c«’'  * 


exP<-naYa1aT,exp(“Va1bl) 


i 
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Since  ya  a  and  ya  b  are  quantities  averaged  by  the  process  described 
in  Equation  (25),  the  l/V  factor  can  be  brought  out  explicitly  to  write 


C(t)  =  exp[-(paYaia  +  pbYaib)Tl 


(48) 


where  pa  and  are  the  molecular  number  densities  of  the  absorber 
and  perturber  gases  respectively.  By  using  this  correlation  function 
in  Equation  (30),  the  interruption  broadened  line  shape  is 


>(Au) )  =  - 


cNLC'“"'~’ 


(49) 


R  I 

where  a  is  the  real  part  of  Ya  a  +  b  and  a  is  the  imaginary 

part  of  the  same  function.  Thus,  the  function  which  describes  the 
near  line  center  profile  is  the  Lorentz  line  shape.  The  term  cJ  repre¬ 
sents  a  line  shift  caused  by  the  impact  approximation  collisions,  while 

r 

a  represents  the  halfwidth  of  the  absorption  time. 

The  derivation  of  Equation  (49)  for  the  long  correlation  times 

implied  by  the  impact  approximation  did  not  explicitly  depend  on  the 

functional  forms  of  the  sel f-broadened  and  foreign-broadened  interaction 

Hamiltonians  written  in  Equations  (37a)  and  (37b).  Instead,  these 

interaction  potentials  influence  the  nature  of  y  and  Y  ■,  respectively 

a^a  a^o 

within  the  derived  line  shape. 

The  pressure  dependence  of  the  line  width  and  line  shift  terms 


in  Equation  (49)  arc*  already  apparent  since  the  molecular  number  densities 
pa  anr*  Pb  are  PorPurt'ona'  t°  pressure  by  the  ideal  gas  law. 
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a  =  AlPa  f  A2pb 
a1  =  BlP,  +  B2pb 


(50) 


where  A^,  A^,  Bj,  arid  B2  are  constants  related  to 
is  the  partial  pressure  of  the  absorber  and  p^  is 


Y,  and  v  ,  ,  p 
Y<V  'axb  Ha 

the  partial  pressure 


of  the  broadener. 


The  temperature  dependences  of  the  line  width  and  line  shift 
are  more  difficult  to  derive.  Here,  the  explicit  forms  for  the  inter¬ 
action  Hamiltonians  in  Equations  (37a)  and  (37b)  become  important. 

A  second  order  perturbation  expansion  of  the  y  and  y  .  [20,21 1 

a  a j  D 

yields  a  temperature  dependence  of  the  form 


(51) 

However,  this  is  not  always  observed  [24,251. 

The  approach  which  has  been  adopted  here  is  the  phase  shift  method 
[26,27],  which  consists  of  only  phase  modulation  of  the  eigenstates, 
but  to  all  orders.  This  approach  leads  to  a  temperature  dependence 
for  the  line  width  of 

a  "y  (52 ) 

for  the  dipole-dipole  interactions  modeling  sel f-broadening 
1 

ot- - 

-j-0.83 

for  dipole-quadrupole  roi  f ore ign-broaoernnq. 


(53) 
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Statistical  Broadening 

Statistical  broadening  is  characterized  by  short  correlation 
functions  caused  by  strong  binary  interactions.  These  interactions 
contribute  to  the  far  wing  absorption  of  a  spectral  line  profile.  No 
perturbation  method  can  be  applied.  However,  the  phase  shift  contri¬ 
butions  are  solvable.  In  addition,  corrections  will  be  included  to 
account  for  energy  level  shifts.  The  time-development  operator  Sca  ^(x) 
can  be  written 


Sca  b^T)  =  exP( ~ 1  n(  t)  ) 


where 


n(x)  =£ 


and 


V*’ 


t  +  T 
r  0 

t  dt  Ha  hW 
tQ  axb 


Vj/(R(t)IJ. 


(54) 

(55) 

(56) 


The  parameter  m.  is  the  power  of  the  exponent  in  the  Hamiltonians  of 
Equations  (38a)  and  (38b).  For  H^O-H^O  interactions  mj=6,  while  for 
H2O-N2  or  H^O-Op,  nOj  =  4.  The  correlation  function  initially  written 
in  Equation  (26)  can  now  be  expressed 


Ca"b(T> 


LS 


Ip  «usb|exp(  i  n(  x)  lusb><stsb  [exp(-i  n(  t)  1 8.s  b  >  >R 
b  sb  o,v,to 

(57) 


This  can  be  written 


nK 


1+pb  .I  p  «usb|exp(in(  x) )  jus^><)tsb  |exp(-i  n(  x))|£sb>-l>R 


sb'  sb 


R0*v*t0 


where  <> 


R  .  v ,  t 
0  ’  0 


r<  >  R  .  and  =  n^/V.  The  fact  that  np  is 

r\  ’  *  L  n 


very  large  permits  the  correlation  function  to  be  written  C28 1 


caib(T)  =  exp 


y  P  h  <1-  <usb  | oxp ( i  n!  t)  )  |usb  ><9sb  |exp(-i  n(  x)  )|ts 


(59) 


•  <b>' 


R„v,t„ 
0  ’  0 


J 
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The  matrix  elements  in  Equation  (59)  can  be  manipulated  to  yield  the 
form 

<usb|  eln|  usbx;tsb|  e_inj£sb>= 

e1*1  u"n£^+einu  ft£  +  en^u  +  ft^  (60) 

where 


nu  =  <usb|  n|  usb>. 


n£  =  <«.  sb|  n|  s.sb>. 


u  L  ^V—Y.  h,  <usb|n|  u'sb  >  ...  <u(n-1)sb(n-D|n| 
*  H'sD 


4(n-l)sb(n-l) 


and 


r  b,<tsblll!1'sb,>...<t(n-l)sb(,l-l)|n|tsb;. 


^ (n-l )sb(n -1 ) 

The  prime  on  the  summation  sign  indicates  that  all  the  indices  can 
not  be  equal.  Theft's  are  the  level  shift  terms  and  have  no  closed 
form  expression.  Let, 


-  ,  ,  ,nu 
f k(t  )  =  e  ft  +  e 


■in. 


aA- 


(61) 
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Now  the  correlation  function  becomes 
n,_ 


Ca,b(x)  =  exP(-PbI'P  J<Ue 

1  sb  s 


l(VV>Ro,v,to  + 

(62) 


i(nu-nJ  ’ 

The  phase  shift  term  <l-e  >  D  .  can  be  solved.  The  major 

contribution  to  the  phase  term  comes  in  the  region  t  <  Rq/V  [29]  where 
v  is  the  velocity  of  the  colliding  molecule.  The  sum  over  t  implied 
by  the  bracket  in  Equation  (62)  as  defined  in  Equation  (25)  is  con¬ 


verted  to  an  integral  which  gives  its  major  contribution  between  -R  /v 
and  Rq/v.  In  this  region 

-bi.i  i  ..  _  b 


/>At 


m --even 

J 


(63) 


where  x  represents  the  subscripts  l  and  u.  The  integral  over  t  gives 


1-(nu-nJ  1  s  . 

<l-e  l>R0,v,t0 


Rn/ v  •  i  \  i 

0  Ha.-nJ 

dt  ( 1-e  u  *)> 

R  R„,v 

0  o’ 


or 


,  J(V\}  '  2Ro  n  i(vVx 

<1_e  >Ro'»  ' 


(64) 


Evaluating  the  real  and  imaginary  parts  separately  we  have 
2R 

<  (l-cos(n  -n  ))>R  = 

V  o’ 


4rr  /t\  3/m-  3/m  -  3/m 

mj  W  (v  ju-  V  *j  r<-  HT  >  C0S  ImT  >  "reven 

J  J 


(65) 
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The  general  correlation  function  from  Equation  (62)  can  now  be  written, 

n,  3/m.. 

C^bU)  =  expC-Ph  l  Psb  {(Yj+ie)i  J  +<fb(  t)  >Ro,v,to  ^  ^  (67) 

Separating  real  and  imaginary  terms  in  the  exponent, 

n.  3/m- 

Cajb(T)  =  Gxp[-pb^  psbV  +  Re[<Fb(^  Ro,v,t0  11X 

r  v  3/m- 

exp[-f|P„  >.  »sbCT  -  l"<VO>RoiV>tol].  ,68) 


When  a  correlation  function  of  this  form  is  substitued  into  Equation 
(32),  the  expression  in  the  second  exponential  combines  with  the  ui 
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to  produce  a  line  shift  term  in  the  final  line  shape  jc(Aw).  It  is 
the  amplitude  of  the  correlation  function  which  yields  the  line  profile 

nb  3/m, 

!calb('r)l  =  exp[-xbjgbj(x)T  J]  (69) 

where  the  substitutions 


=  *  l  °b  TJ 

s°  S 


95j(t) 


=  1  + 


Re  [<f h(  t)  >r  tl 
_ b _ V^o_ 

3/ in  - 
YjT  J 


(70a) 


(70b) 


were  made. 

The  general  correlation  function  for  a  binary  mixture  is  then 


3/m.  3/m, 

Cfw(t)  =  expf-^.  gbj(T)T  JJexp[-Aakgak(T)T 

and  the  far  wing  line  shape  is 


(71) 


JcFW^^  =  1  ]  dx  exp  [-i  Au)t]Cfw(  |)| 


(72) 


An  asymptotic  expansion  of  the  exponentials  in  Equation  (71)  yields 

3/m, 


"k  c  3/m,  1 

c  ,T)  .  y  y  ('xak9ak^^  >s  <~Wt)t  J) 

°PVJ\  '/  {‘  r  •  ~  t  r 

s=0  t-0 


(73) 


This  yields  a  far  wing  line  shape  which  is  written  as 
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JCFw(Aaj)  TT  Re^ 


CO  CO 

l  l 

s=0  t=0 


3/m  ■ 


J] 


Jo 


3/mk  S 

dr  [exp  -iZxajt  ][-gak(x  )t  ] 


(74) 


Making  the  substitution  0=Awt,  the  far  wing  line  shape  becomes 


J'CFW^) 


oo 


j  Am| 


1 


+ 


3t 


x 


Re 


co 


d0  exp(-iO) 
.  Jo 


(75) 


The  leading  terms  in  this  expression  occur  for  s=l,  t=0  and  for  s=0, 
t=l  since  the  single  term  where  both  s  and  t  are  zero  is  zero.  These 
two  leading  terms  have  the  physical  interpretation  of  being  the  self- 
broadened  and  foreign-broadened  contributions  to  the  far  wing  line 
shape.  If  we  retain  only  these  leading  terms. 


W4“)  ="  uV37”/  Rei< 


*  Jul- 


M1,3/mj 


Re 


d0exp(-iO ) 

ct)  exp(-io) 


‘  gak(/Vo 


bjy 


0- 

Aio 


3/m . 


3/m, 


1} 


(76) 


Evaluation  of  the  integrals  requires  substitution  of  Equations  (61) 

and  (70b).  In  general,  it  can  be  argued  that  the  resulting  functions 

will  be  different  when  m,-m.. 

*  J 
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The  far  wing  line  shape  can  be  written 


=  - 


^ak 
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1+3/m,  "k 


-  T 


+  XbJ 


Tm  K,Mj) 


Mltvrai  V  * 


(77) 


where  =6  for  self-broadening  and  m ^ =4  for  foreign-broadening  by 
or  02. 

The  statistical  line  shape  represented  by  Equation  (77)  is  similar 
to  the  results  presented  by  Fomin  and  Tvorogov  [30].  In  their  development, 
the  term 


1 


IH 


l+3/m 


is  multiplied  by  an  exponential  factor  which  is  proportional  to 


m-1 

{exp(-j  Aw|  m  )2/3}. 


For  self-broadening,  this  factor  in  their  line  shape  becomes  exp(|  Auj|  ^  ) > 
while  for  foreign-broadening  this  factor  is  exp([  Awj  <3'3<3)  •  These  fac¬ 
tors  are  expected  to  be  present  in  the  T-functions  of  Equation  (77). 

The  p-function 

In  the  previous  two  sections,  interruption  broadening  and  statistical 
broadening  were  presented.  The  total  line  shape  is  constructed  by 
combining  these  two  resulting  line  shapes.  We  have  chosen  to  do  this 
by  the  application  of  band-pass  and  band-stop  filter  functions  to  the 
two  derived  line  shapes.  The  total  line  shape  is  written 

Jc(Av)  =  N  jCNLC(A.)p(v)+jCFW(Av)(l-p(s;)) 


(73) 
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where  we  have  made  a  change  of  variable  from  frequency  w  (measured  in  radian/ 
sec)  to  frequency  v  (measured  in  cm'").  The  normalization  constant 
N  assures  that  Equation  (2)  is  satisfied.  The  function  p(y)  is  a  nor¬ 
malized  band-pass  filter  function  centered  at  vQ.  Guidelines  for  choosing 
a  suitable  p-function  are  given  in  Table  2  . 

Based  on  these  guidelines,  the  p-function  was  chosen  to  be 


P(v) 


1 

2 


0 


K 

which  is  shown 


in  Figure  6. 


j  Avj  <  5  cm 
|  Av|  5  cm' 


(79) 


I 

l 

Ficure  6.  p(  v)  versus  Au 
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TABLE  2 


1.  The  p-function  must  be  equal  to  1  at  v  =  vQ. 

2.  The  p-function  must  be  zero  for  all  v  greater  than  a  cut-off  fre¬ 

quency  v,. 

a 

3.  The  p-function  must  be  a  monotonically  decreasing  function  in  the 
region  vQ  <  v  <  \>a. 

4.  The  p-function  must  be  symmetric  about  y  . 

5.  The  derivative  of  the  p-function  must  be  zero  at  y  and  at  v  =  v  . 

6.  It  is  desirable  for  the  total  line  shape  as  expressed  by  Equation 
(78)  to  be  easily  integrated  to  evaluate  the  normalization  constant 
N. 


The  cut-off  at  5  cm-1  from  the  line  center  was  chosen  so  that 
the  p-function  would  remain  essentially  unity  within  several  halfwidths 
of  the  water  vapor  line  center.  This  allows  interruption  broadening 
to  dominate  near  line  center.  The  stop-band  filter  function  (l-p(v)) 
allows  the  stati-ticcl  line  shape  to  dominate  beyond  5  crrf^  (50  to  10G 
halfwidths  from  lin  center). 

The  region  whe:  e  p(v)  and  l-p(v)  are  approximately  equivalent 
is  the  intermediate  wing  region.  This  method  permits  a  smooth  blending 
■  interruption  broadened  and  statistical  broadened  line  shapes, 

’  ’•iif  >•,  the  necessity  for  solving  or  approximating  the  inter- 

<*  1  i  ni  ten. 

■  •  •  a  the  p-function  and  the  exact  choice  of  cut- 
important  to  the  study  of  the  water  vapor 
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continuum  problem  in  the  infrared.  The  reason  for  this  is  that  the 
continuum  absorption  is  theorized  here  to  be  caused  by  the  overlap 
of  the  far  wings  of  the  strong  water  vapor  absorption  lines.  The 
intermediate  wing  region  and  near-1 ine-center  region  enter  the  problem 
only  indirectly  in  that  they  must  be  known  to  normalize  the  total  line 
shape  function  properly.  Furthermore,  the  p-function  formalism  permits 
the  blending  of  the  different  pressure  and  temperature  dependences 
of  the  interruption  broadened  and  statistical  broadened  line  shapes 
as  will  be  seen  later  [15]. 

CONCLUDING  REMARKS 

A  total  line  shape,  valid  for  water  vapor  absorption  lines  through¬ 
out  the  infrared  has  been  derived  by  proper  consideration  of  both  impact 
and  statistical  collisional  processes.  Self-broadening  is  assumed 
to  occur  through  the  electric  dipole-dipole  interaction,  while  foreign- 
broadening  is  assumed  to  occur  through  the  electric  dipole-quadrupole 
interaction.  Thus,  the  derived  shape  is  suited  for  atmospheric  (N^ 
or  0^)  broadening  of  H^O  absorption  lines. 

Several  approximations  were  made  in  the  process  of  deriving  a 
proper  infrared  line  shape  for  H?0  molecules.  Fundamentally,  the  binary 
collision  approximation  and  the  adiabatic  approximation  for  the  energy 
levels  in  the  absorber  molecule  restrict  the  temperature  and  pressure 
regimes  over  which  this  spectral  line  shape  might  be  valid. 

In  the  next  section,  implementat  on  of  the  total  hm  shape  will 
be  discussed,  and  the  ability  of  this  line  shape  to  predict  continuum 
absorption  in  the  8-12  mn  and  in  the  3-4  pm  atmospheric  transmi ttance 
windows  will  be  discussed. 
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Part  II:  Implementation  of  the  Line  Shape 


INTRODUCTION 

In  the  preceding  paper  we  developed  an  expression  for  the  pressure- 
broadened  absorption  line  shape  of  l^O  molecules.  Only  binary  collisions 
were  assumed,  and  both  self -broadening  and  foreign-broadening  with  were 
included.  The  spectral  line  shape  was  developed  by  first  considering 
the  time  evolution  of  the  correlation  function  of  the  energy  levels  of 
the  colliding  molecules.  Collisions  were  categorized  by  their  severity: 
weak  collisions  contributed  only  to  the  near-line-center  of  the  spectral 
profile,  while  strong  collisions  contributed  to  the  line  shape  everywhere 
including  the  far  wings.  The  study  of  H^O  continuum  absorption  must  in¬ 
clude  a  study  of  the  accumulation  of  these  far  wings. 

Statistical  broadening  [ l3  was  used  in  our  development  to  formulate 
the  far  wing  component  to  the  spectral  line  shape.  Recently,  another 
line  shape  development  used  only  the  impact  approximation,  modified  to 
account  for  variations  in  collision  duration  t.2 J .  It  is  not  clear  at 
this  time  whether  there  is  any  significant  difference  in  the  line  shapes 
produced  by  these  two  separate  approaches. 

In  this  paper  we  demonstrate  how  the  spectral  line  profile,  defined 
by  Equations  (49),  (77),  and  (78)  in  the  preceding  paper,  is  implemented. 
Furthermore,  comparison  between  experimental  absorption  data  and  the  ab¬ 
sorption  predicted  by  this  line  shape  is  given.  We  will  show  that  a  line 
shape  formulation  is  adequate  to  predict  all  of  the  behaviors  of  water 
vapor  continuum  absorption  over  the  ranges  of  pressure  and  temperature 
of  atmospheric  interest.  Data  will  be  shown  for  the  8-14  urn  and  3.5-4 
urn  atmospheric  transmission  windows. 
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HISTORICAL  PERSPECTIVE 

Figure  1  demonstrates  the  inability  of  the  Lorentz  line  shape  to 
model  the  observed  pressure  dependence  of  water  vapor  continuum  absorption 
near  10  urn  for  H^O  samples  pressure  broadened  to  1  atm  with  N^.  The  data 
were  recorded  using  a  frequency-stabi 1 ized  CC^  laser  and  a  long-path  ab¬ 
sorption  cell  described  in  several  references  [35].  The  curve  marked 
A  was  generated  from  data  on  the  AFGL  tape  C6 1  using  a  bound  of  20  cm*'1. 
That  is,  for  the  calculations  at  each  H2O  partial  pressure,  only  the 
H2O  absorption  lines  within  a  range  of  -20  cm~^  from  the  measurement  fre¬ 
quency  were  assumed  to  contribute  to  the  absorption  at  that  frequency. 

All  other  H^O  absorption  lines  were  assumed  to  be  too  far  away  to  make 
a  significant  contribution.  This  bound  of  20  cm-1  became  a  standard  in 
early  line-by-line  calculation  programs  because  it  was  believed  to  be 
a  reasonable  compromise  between  spectral  representation  and  computing 
time. 

The  curve  marked  B  in  Figure  1  was  generated  by  accumulating  the 
contributions  of  al 1  H^O  absorption  lines  from  0  to  5000  cm'  .  The  con¬ 
tributions  at  the  measurement  frequency  were  computed  using  the  Lorentz 
line  shape  for  all  lines  listed  on  the  data  tape.  At  low  H^O  partial 
oressures,  this  line  shape  overestimates  the  observed  absorption,  while 
at  high  H^O  partial  pressures  it  underestimates  the  observed  absorption. 
Despite  this  failing,  the  difference  between  curve  A  and  curve  B  demon¬ 
strates  that  significant  contribution  to  the  total  absorption  can  be  found 
in  the  accumulation  of  the  far  wings  of  strong  absorption  lines.  The 
Lorentz  profile,  however,  is  not  an  adequate  model  for  the  far  wings  of 
lines,  so  it  is  not  surprising  that  curve  B  fails  to  provide  suitable 
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Figure  1.  Pressure  dependence  of  water  vapor  continuum  absorption. 
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agreement  with  the  measured  absorption.  Motivated  by  the  hope  that  an 
accurately  derived  HpO  absorption  line  shape  would  improve  the  modeling 
of  continuum  absorption,  we  undertook  the  work  presented  in  this  paper 
and  in  the  previous  section. 

THE  LINE  PROFILE 

As  mentioned  in  the  last  section,  the  absorption  coefficient  at  fre¬ 
quency  v  (in  cm-'*')  can  be  written 

k(v)  =  I  Si  j(v-vo_;8i)  (1) 

where  is  the  strength  of  the  ith  absorption  line  and  j(v-v  ;8-)  is 
the  shape  of  the  absorption  line  centered  at  and  which  has  a  set  of 
shape  parameters  { 8-}  .  The  major  features  of  the  total  line  shape  have 
been  derived  in  the  previous  section.  It  was  found  that 


jc(Av)  =  N{jCNLC(  Av)p(  Av)  +  JCFW(  Av)(l-p(  Av))  }  (2) 


where  the  subscript  C  denotes  that  the  line  shape  is  collision  broadened, 
and 


R, 
a  2 


Jcnlc<  4v)  =  -  (WiQf 


(3) 


%W(Av)  =  “ 


ak 


T  _  (in.  ,Av  ) 


x  hi 

+  — “--1.  0/  ‘  T  m  ( m  - ,  Av) 

1+3/m  •  m j’ 

| Av|  J  J 
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1/2  +  1/2  cos 
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p(Av)  = 


Av) 


| Av|  <5  cm 


-1 


|  Av|  2.  5 


cm 


-1 


(5) 


N  is  the  normalization  constant.  Normalization  requires  that 


fj  (Av)dv  =  1.  (6) 

Jo 


In  order  to  eliminate  the  singularities  in  the  far  wing  line  shape  func¬ 
tions  at  Av=0,  we  now  adopt  the  Varanasi  formulation  [7]  which  includes 
an  additive  term  in  the  denominators  of  Equation  (4).  Furthermore,  we 
now  set  m(<=6  and  mj=4  as  discussed  in  the  previous  paper. 

r 


1 


0.3198X 


aJL 

3 

Aa6 


T6H=6’Av) 


i 


+ 


0.4334Ab4 
I Av| 1 *75+( A54) 


7/3 


1 

T4(mj=4,Av) 


(7) 


This  modification  of  the  line  shape  function  does  not  significantly  alter 
the  value  of  the  function  when  Av  is  very  large.  It  only  assures  that 
this  function  is  well-behaved  at  line  center  for  normalization  purposes. 
Furthermore,  in  our  formulation  the  p-function  is  turning  off  the  far 
wing  line  shape  near  Av=0,  so  any  constant  value  of  at  Av=0  arti- 

fically  introduced  by  this  modification  is  given  a  weighting  of  zero. 

The  multiplicative  constants  which  appear  in  Equation  (7)  are  aids 
in  normal ization  introduced  by  the  Varanasi  formulation.  An  evaluation 
of  ^a6  an(^  \)4  'n  ^ar  w’n9  data  (Av>>^a6’*b4^  wou^  simply  be  scaled  by 
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these  constants. 

Now,  the  functions  Tg  (mk=6,  Av)  and  T4  (rrtj=4,Av)  must  be  discussed. 

As  mentioned  in  the  previous  section,  Fomin  and  Tvorogov  [8j  derived  a  term 

of  the  form 


"  m-r 

2/3  \ 

G  1  Av|m 

T  j 

J 

where  G  is  a  constant  and  T  is  the  temperature.  An  exponential  decay 
of  some  form  is  expected  in  the  far  wings  of  absorption  lines  as  discussed 
by  Clough  and  coworkers  [2].  For  self-broadening  (m=6)  the  above  term 
is  proportional  to  exp(-|Av|u*  ),  and  for  foreign-broadening  this  term 
is  proportional  to  exp( - | Av|^*50) . 

We  have  adopted  a  simplification  of  this  term.  We  have  set 

T6(mk=6,Av)  =  exp(-Ga6|Av|0'5)  (8a) 

and 

T4(mj=4,Av)  =  exp(-Gb4|Au|0’5).  (8b) 

We  have  given  the  same  frequency  dependence  to  both  the  self-broadened 
and  foreign-broadened  terms  in  order  to  facilitate  the  normalization  pro¬ 
cess.  The  origin  of  these  terms  is  the  level  shift  contribution  which 
was  not  fully  incorporated  in  our  phase  shift  approach.  Such  terms  are 
difficult  to  derive  by  perturbation  methods  since,  by  definition,  in  the 
far  wings  of  infrared  absorption  lines  the  perturbation  energy  (Av)  is 
no  longer  small  compared  to  the  unperturbed  energy  levels  of  the  absorbing 
molecule. 

The  G  parameters  above  are  written 


r  (  V 296^  Ya(V^ 

Ga6  =  a9(v)[^-TrJ 
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(9a) 


and 


where  the  parameters  9a(v),  gb(v),  Ya(v),  and  ^^(v)  depend  on  the  vibra¬ 
tional  quantum  numbers.  In  particular,  it  can  be  anticipated  that  the 
parameters  might  be  slightly  different  for  the  pure  rotational  band  (v=0) 
where  various  rotational  levels  are  more  easily  accessible  through  thermal 
energy  than  in  the  case  for  the  vibrational  energy  levels.  Thus  we  might 
expect  one  set  of  parameters  to  describe  the  shape  of  lines  in  the  ro¬ 
tational  band  and  another  set  of  parameters  to  describe  the  shape  of  ab¬ 
sorption  lines  of  the  vibrational  bands. 

Equation  (7)  is  then  written 


(ID 
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where  the  terms  have  been  defined  previously.  For  self-broadening,  this 
can  be  written 


^  ^  ^  >  (v6u(T)-v„(T))1/2  t‘/2r(- cos  (  J) 

(12) 

The  potential  for  self-broadening  is  seen  in  Equation  (40)  of  the  previous 
paper  to  be  proportional  to  1/T.  Therefore  we  can  write 


where  pg  is  the  absorber  pressure  and  Al(v)  is  a  parameter  which  should 
have  some  dependence  on  the  vibrational  quantum  numbers,  and  296  is  the 
reference  temperature. 

The  foreign-broadening  parameter  can  be  written 
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which  can  be  written  simply  as 

Xb4  *  A2(v,J)(^Pb  (15) 

where  p^  is  the  broadener  pressure  and  A2(v,J)  is  a  parameter  which  de¬ 
pends  on  the  vibrational  quantum  number  and  on  the  rotational  quantum 
number  J.  8ecause  of  the  statistical  averaging  over  rotations,  the  HpO- 
HpO  potential  matrix  elements  for  strong  collisions  are  independent  of 
angle  and  not  a  function  of  the  rotation  quantum  number  J.  For  Np-HpQ 
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however,  the  interaction  times  are  shorter  than  the  rotational  times; 
thus,  no  averaging  over  rotations  is  performed.  Thus,  A2(v,J)  will  be 
a  function  of  the  J  quantum  number.  This  J  dependence  can  be  expressed 


in  terms  of  the  J  dependence  of  the  near-line-center  halfwidth  aQ  [9] 


A2( v, J)  =  A2( 


^(or) 


1.125 


where  aQ  is  the  near-1 ine-center  halfwidth,  and  0.07  represents  a  mean 
value  of  aQ  for  pressure  broadened  water  vapor  lines  in  the  mid-infrared. 
The  pressure  and  temperature  dependence  of  the  halfwidth  a  was  taken 


where  pa  is  the  absorber  pressure  and  pa  is  the  broadener  pressure.  B 

is  the  self-to-foreign  broadening  ratio  which  is  typically  set  to  5. 

The  self-broadened  component  of  the  halfwidth  has  a  temperature  dependence 

of  1/T,  while  the  foreign-broadened  component  has  a  temperature  dependence 
0  83 

of  1/T  '  as  discussed  in  the  previous  section. 

In  order  to  write  the  absorption  coefficient,  the  line  strength  is 
needed. 


The  term  written  as  Sgg^  is  a  correction  to  the  line  strength  term  which 
takes  account  of  the  fact  that  for  shifted  energy  levels  the  Boltzmann 
population  distribution  is  not  calculated  from  the  energy  level  separation 
hvQ,  but  rather  from  the  total  energy  separation  hvQ+AE.  The  correction 
is  important  for  the  rotational  band  where  different  rotational  levels 
are  populated  at  room  temperature.  The  correction  is  far  less  important 


3« 


for  the  vibrational  bands.  The  term  SCQR  can  be  shown  [9]  to  take  the 
form 


5C0R 


(l-exp(-hcv/kT ) ) 

tt  Ti - r: — near  line  center 
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Since  these  terms  are  frequency  dependent,  we  can  include  them  in  the 


line  shape  term.  These  terms  only  affect  absorption  lines  below  approxi¬ 
mately  1000  crn~ ^ .  Now  the  total  line  shape  can  be  written 
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As  mentioned  earlier,  the  near-line-center  profile  and  the  far  wing 
profile  were  blended  through  the  intermediate  wing  region  by  the  use  of 
the  p-function. 

k(v)  =  I  Si  ( JcNLC^ v_vo  )p(Av)  + 
i  i 

%w(v-vo  ;Bi)(l-p(Av))}  (14) 

i 

The  normalization  constant  for  each  line  is  determined  by  imposing  Equa¬ 
tion  (6)  on  each  line. 

COMPARISON  WITH  EXPERIMENT 

A  computer  program  was  written  which  incorporated  this  line  shape 
in  a  line-by-line  calculation  of  absorption  in  the  infrared. 

Before  calculations  of  the  contribution  of  far  wings  to  H^O  continuum 
absorption  could  be  made,  the  various  parameters  had  to  be  evaluated. 

This  was  done  by  fitting  predicted  absorption  coefficients  to  measured 
absorption  coefficients  in  spectral  regions  where  line  absorption  clearly 
dominates.  It  must  be  mentioned  here  that  absorption  data  in  the  window 
regions  where  continuum  can  be  measured  was  not  used  in  this  evaluation 
of  the  parameters.  The  regions  chosen  for  the  parameter  evaluation  were 
located  in  the  absorption  bands  where  absorption  coefficient  on  the  order 
of  100  km~^  are  typical.  Since  continuum  absorption  is  on  the  order  of 
0.1  km~*  in  the  window  regions,  any  possible  residual  continuum  absorption 
in  the  regions  where  the  line  parameters  were  evaluated  would  have  been 
a  small  fraction  of  the  total  absorption. 

Data  provided  by  Burch  [10]  were  used  to  evaluate  the  parameters 
which  define  the  rotational  band.  Data  recorded  with  a  CO  laser  at  our 
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laboratory  [9]  were  used  to  evaluate  parameters  of  the  v  vibrational 
band  near  6.6  wn,  and  HF  laser  measurements  made  by  White  et  al  [11]  were 
used  to  estimate  the  parameters  in  the  spectral  region  near  3  urn. 

In  all  cases  it  was  clear  that  insufficient  temperature  data  existed  to 
make  an  accurate  estimate  of  certain  temperature  controlling  parameters 
such  as  Ya(v)  and  Y^v).  For  these  parameters  it  was  decided  to  vary 
them  only  slightly  if  necessary  from  the  value  0.67  given  by  Fomin  and 
Tvorogov  [8] . 

No  attempt  at  a  least  square  analysis  was  made  since  the  data  do 
not  yet  possess  the  accuracy  and  completeness  to  make  such  an  analysis 
meaningful.  Table  1  lists  the  resulting  parameters  for  each  H^O  absorp¬ 
tion  band.  As  can  be  seen,  the  largest  deviation  between  parameter  values 
occurs  between  the  rotational  band  and  the  rest  of  the  absorption  bands. 
Although  this  trend  was  anticipated,  no  physical  significance  should  be 
placed  on  the  small  differences  which  appear  to  exist  between  the  paramet¬ 
er  values  for  these  bands.  Again  it  is  important  to  emphasize  that  absorp 
tion  data  in  the  window  regions  were  not  used  in  the  evaluation  of  these 
parameters. 

Once  the  parameters  were  estimated,  we  used  the  total  line  shape 
to  calculate  the  absorption  at  CO,,  laser  frequencies  in  the  10  urn  atmos¬ 
pheric  transmission  window.  Figure  2  shows  the  results  of  the  calculation 
at  room  temperature  for  the  pressure  dependence  of  the  absorption  coef¬ 
ficient  at  one  of  the  laser  lines.  A  bound  of  1000  cm~*  was  used,  al¬ 
though  we  found  that  increasing  this  made  no  significant  difference  to 
the  calculated  absorption  coefficient.  The  predicted  absorption  is  in 
very  good  agreement  with  the  measured  absorption  at  all  pressures.  It 
should  be  pointed  out  that  the  curvature  in  the  data  (seen  as  a  p  depend- 
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Figure  2.  Room  temperature  modelling  of  water  vapor  pressure  depende 
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ence)  is  well  modelled  by  this  line  shape  theory.  This  is  a  manifestation 
of  the  increased  importance  of  H^O-F^O  interactions  over  F^O-^  interac¬ 
tions  for  strong  collisions  as  developed  in  the  previous  paper.  Using 
the  Lorentz  profile  as  in  Figure  1  the  curvature  in  the  data  was  not  well 
model  led. 

Our  results  in  modelling  the  overall  continuum  absorption  throughout 
the  8-14  un  atmospheric  window  can  be  seen  in  Figure  3.  Measurements 
of  the  self-broadened  coefficient  C$  made  by  Burch  [12]  were  compared 
to  predictions  calculated  with  various  classic  line  shape  formalisms. 

Using  our  line  shape  formalism,  a  very  good  fit  to  the  data  was  achieved. 

We  have  also  applied  the  line  shape  theory  to  the  3.5-4  pm  trans¬ 
mission  window.  DF  laser  measurements  made  by  White  et  al  [13]  did  not 
confirm  earlier  spectroscopic  measurements  made  by  Burch  [14].  White's 
measurements  were  consistently  higher  than  the  absorption  values  reported 
by  Burch,  and  the  curvature  as  a  function  of  frequency  in  White's  data 
throughout  this  region  was  much  greater  than  seen  by  Burch. 

Figure  4  shows  the  results  of  our  calculations  in  the  3.5-4  pm  region. 
The  squares  are  the  total  absorption  calculated  using  water  vapor  absorp¬ 
tion  lines  and  a  partial  pressure  of  15  torr  in  an  environment  of  745 
torr  ^  This  spectral  region  contains  many  absorption  lines  of  the  HDO 
molecule.  Burch  was  careful  to  avoid  regions  close  to  the  centers  of 
these  absorption  lines  when  he  made  his  measurements.  We  can  accomplish 
this  in  our  calculations  by  subtracting  the  local  line  contribution  within 
a  few  wavenumbers.  The  triangles  show  the  residual  continuum  absorption 
predicted  by  our  line  shape  calculation.  The  agreement  with  Burch's  data 
is  quite  good. 
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Figure  3.  Experimental  results  in  the  8  to  14  pm  region. 
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Modelling  the  temperature  dependence  is  much  more  difficult  because 
a  precise  understanding  of  intermolecular  interactions  does  not  exist 
(particularly  for  the  rotational  band)  and  accurate  temperature  data  do 
not  exist  for  all  of  the  absorption  bands  of  H^O  in  the  infrared.  With 
the  values  of  the  line  shape  parameters  listed  in  Table  I,  we  calculated 
the  absorption  coefficient  at  the  P(20)  CO2  laser  line  at  944.1945  cm’^. 
Data  which  exist  [3,15]  on  the  absorption  in  this  region  indicate  a  strong, 
negative  temperature  dependence.  For  this  study,  values  of  Ya  and 
were  fixed  to  0.67.  Improved  data  in  all  spectral  regions  is  necessary 
before  refinement  of  these  parameters  can  be  made. 

Figure  5  shows  the  results  of  our  modelling  of  the  temperature  de¬ 
pendence.  A  negative  temperature  dependence  is  predicted  from  our  line 
shape  model.  Previously,  the  negative  temperature  dependence  of  the  ab¬ 
sorption  was  taken  as  strong  evidence  of  dimer  absorption. 

An  adequate  line  shape  model  for  atmospheric  H^O  is  important  not 
only  for  modelling  the  continuum  absorption,  but  also  for  calculating 
the  effects  of  overlap  between  F^Q  absorption  line  features  and  weak  ab¬ 
sorptions  from  a  possible  trace  gas  in  the  atmosphere.  Recently,  Menyuk 
et  al.  [16]  used  the  differential  absorption  lidar  technique  to  detect 
trace  amounts  of  NO  in  the  atmosphere.  Interference  from  strong  absorp¬ 
tion  lines  of  H^O  were  not  adequately  compensated  for  by  two  separate 
programs  provided  by  AFGL.  We  were  asked  to  calculate  the  absorption 
coefficient  of  pressure-broadened  water  vapor  at  their  measurement  frequen¬ 
cies.  Using  the  results  which  we  provided  to  them,  they  recomputed  the 
effects  of  H^O  interference  in  their  experiment  and  found  much  better 
agreement  with  their  experimental  results. 
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Figure  5.  Comparison  of  experimental  and  theoretical  temperature  depend¬ 
ence  of  Cs  dt  944.1945  cm'1. 
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CONCLUSIONS 

In  these  two  papers  we  have  developed  and  discussed  the  shape  of 
pressure-broadened  H^O  absorption  lines  in  the  infrared  spectral  region. 
The  water  vapor  continuum  in  both  the  8-14  um  and  in  the  3.5-4  pm  atmos¬ 
pheric  windows  is  quite  well  modelled  by  this  absorption  line  shape. 

The  success  of  the  modelling  extends  to  a  prediction  of  the  negative  temp¬ 
erature  dependence  observed  near  10  Pm. 


100 


REFERENCES 


1  R.G.  Breene,  The  Shift  and  Shape  of  Spectral  Lines,  Perqamon  Press, 
1961. 

2  S.A.  Clough,  F.X.  Kneizys,  R.  Davies,  R.  Gamache,  and  R.  Tipping, 
"Theoretical  Line  Shape  for  H^O  Vapor:  Application  to  the  Continuum" 
in  Atmospheric  Water  Vapor  edited  by  Adarsh  Deepak,  T.D.  Wilkerson, 
and  L.H.  Ruhnke,  Academic  Press,  1980. 

3  J.C.  Peterson,  "A  Study  of  Water  Vapor  Absorption  at  CO^  Laser  Fre¬ 
quencies  Using  a  Differential  Spectrophcna  and  White  Ceil,"  Ph.D. 
Dissertation,  The  Ohio  State  University,  June  1978. 

4  R.J.  Nordstrom,  M.E.  Thomas,  J.C.  Peterson,  E.K.  Damon,  and  n.K. 
Long,  Appl .  Opt.  17,  2724  (1978). 

5  J.C.  Peterson,  M.E.  Thomas,  R.J.  Nordstrom,  E.K.  Damon,  and  R.K. 
Long,  Appl.  Opt.  18,  834  (1979). 

6  McClatchey,  R.A.,  Benedict,  W.S.,  Clough,  S.A.,  Burch,  D.E.,  Calfee, 
R.F.,  Fox,  K.,  Rothman,  L.S.,  and  Garing,  J.S.,  U.S.  Air  Force  Re¬ 
search  Laborator ies,  AFCRL-TR-73-0096,  Bedford,  Massachusetts  (1973) 

7  Varanasi,  P.,  Chou,  S.,  and  Penner,  S.S.,  J.  Quant.  Spectroc.  Radiat 
Transfer,  8,  1537  (1968). 

8  Fomin,  V.V.  and  Tvorogov,  S.D.,  Appl.  Opt.  12,  584  (1973). 

9  Thomas,  M.E.,  "Tropospheric  Water  Vapor  Absorption  in  the  Infrared 
Window  Regions,"  Report  784701-5,  August  1979,  The  Ohio  State 
University  ElectroScience  Laboratory,  Department  of  Electrical  En¬ 
gineering;  prepared  under  Contract  DAAG-29-77-C-0010  for  U.S.  Army 
Research  Office. 

10  D.A.  Gryvnak  and  D.E.  Burch,  "Infrared  Absorption  by  C0-,  and  FU0," 
Aeronutronic  Publication  U-6417,  Air  Force  Geophysics  Laborator ies. 
Contract  F19628-76-C-0302,  (May  1978). 

11  W.R.  Watkins,  R.L.  Spellicy,  K.O.  White,  B.Z.  Sojka  and  L.R.  Bower, 
Appl.  Opt.  18,  1582  (1979). 

12  Burch,  D.E.,  Aeronutronic  Publication  No.  U-4784,  Semi-Annual 
Technical  Report,  AFCRL  Contract  No.  F19628-69-C-0263,  U.S.  Air 
Force  (1970). 

13  White,  K.O.,  Watkins,  W.R.,  Bruce,  C.W.,  Meridth,  R.E.,  and  Smith, 
F.F.,  Appl.  Opt.  17,  2711,  (1978). 


D.E.  Burch,  D.A.  Gruyvnak,  and  J.D.  Pembrook,  Aeronutronic  Publica¬ 
tion  No.  U-4897,  "Investigation  of  the  Absorption  of  Infrared  Radi¬ 
ation  by  Atmospheric  Gases:  Water,  Nitrogen,  Nitrous  Oxide,"  Air 
Force  Cambridge  Research  Laboratories,  Contract  No.  F 19628-69- 
C-0263,  January  1971. 

V.N.  Aref'ev  and  V.I.  Dianov-K lokov.  Opt.  Spectrosc.  42,  488  (1977). 

N.  Menyuk,  D.K.  Killinger,  and  W.E.  DeFeo,  Appl.  Opt.  19,  3282 
(1980). 


10? 


SECTION  IV 

FOURIER  TRANSFORM  SPECTRA  OF  H20  SAMPLES 

We  have  collected  two  sets  of  high-resolution  absorption  spectra 
of  H^O  samples.  The  first  set  consists  of  four  spectra  of  pure  T^O  sam¬ 
ples.  These  spectra  were  recorded  at  0.04  cm’*  resolution  using 
an  absorption  path  length  of  821  meters  in  a  temperature-control lable 
White  cell  f.auippeo  with  a  double-pass  optical  scheme.  The  four  spectra 
were  recorded  at  the  different  temperatures  and  T^O  pressures  indicated 
in  Table  1.  For  all  four  of  these  spectra  the  spectral  regions  near 
absorption  band  centers  are  totally  opaque.  The  purpose  for  collecting 
the  spectra  was  the  observation  of  H.,0  absorption  line  positions  in 
the  atmospheric  transmi ttance  window  regions.  The  spectral  coverage  of 
these  Fourier  transform  spectra  extends  from  850  cm'*  to  5000  cm'*.  The 
natural  D/H  ratio  exists  in  the  H^O  samples  used. 

The  second  set  of  data  consists  of  four  spectra  which  demonstrate 
the  pressure-broadening  of  a  single  water  vapor  amount  by  nitrogen.  These 
spectra  were  all  recorded  at  room  temperature.  Varying  amounts  of  N^ 
were  added  to  the  H90  sample.  The  resolution  of  the  spectra  is  0.04 
cm’*  even  though  the  pressure-broadened  lines  are  much  wider  than  this 
resolution.  Table  ?  summarizes  the  experimental  conditions.  The  path 
length  for  these  spectra  was  also  821  meters. 

An  atlas  of  pure  H.,0  absorption  lines  and  an  atlas  of  pressure- 
broadened  H^O  absorption  lines  are  being  prepared.  These  compilations 
will  contain  the  H.,0  absorption  lines  which  lie  in  the  atmospheric  trans¬ 
mittance  windows.  Figures  la  and  lb  show  a  typical  page  of  the  atlas  and 
the  line  directory  table,  respecti vely.  The  experimental  conditions  are 


TABLE  2 

EXPERIMENTAL  CONDITIONS  FOR  NITROGEN  BROADENED  H?0  SPECTRA 
(PATH  LENGTH  =  821  METERS) 
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Figure  1.  b.  HgCdTe 
Pure  H 
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given  in  Table  1  with  the  5  Torr  spectrum  at  the  top  of  the  page.  Fig¬ 
ures  2a  and  2b  show  a  sample  from  the  2250-2260  cm'1  region,  as  well  as 
pure  F^O  samples  as  described  in  Table  1.  Figures  3a,  3b,  4a,  and  4b 
show  the  same  regions  from  the  pressure-broadened  data  set  as  described 
in  Table  2.  The  region  from  850  to  5000  cm"1  is  covered  in  each  of  these 
compilations  at  10  cm-1  per  page.  Regions  near  absorption  band  centers 
where  the  large  density  of  strong  lines  causes  total  absorption  are  not 
plotted.  Nevertheless,  each  atlas  of  absorption  spectra  contains  over 
300  pages.  Thus,  it  is  not  feasible  to  include  such  compilations  in  this 
report.  Instead,  we  simply  report  the  availability  of  these  data  sets. 

An  accurate  line  finding  routine  has  been  employed  to  locate  and 
catalog  most  of  the  absorption  lines  in  the  spectra.  Each  atlas  is 
accompanied  by  a  line  listing  which  gives  the  position  (in  cm"1)  of  each 
absorption  line  located  and  numbered  by  the  line  finding  routine.  The 
frequency  accuracy  is  estimated  to  be  about  +  0.003  cm”1. 


Figure  2.  b.  InSb 
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